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Type la Supernova Explosions in Binary Systems: 
The Impact on the Secondary Star and its Consequences 
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ABSTRACT 

One method of discriminating between the many Type la progenitor scenarios is by 
searching for contaminating hydrogen and helium stripped from the companion star. 
However, this requires understanding the effect of the impact of the supernova shell on 
different companion stars to predict the amount of mass stripped and its distribution 
in velocity and solid angle for the types of binary scenarios that have been proposed as 
Type la progenitor models. 

We present several high-resolution 2-D numerical simulations of the impact of a 
Type la supernova explosion with hydrogen-rich main sequence, subgiant, and red 
giant companions. The binary parameters were chosen to represent several classes of 
single-degenerate Type la progenitor models that have been suggested in the literature. 
We use realistic stellar models and supernova debris profiles to represent each binary 
system. For each scenario, we explore the hydrodynamics of the supernova-secondary 
interaction, calculate the amount of stellar material stripped from the secondary 
and the kick delivered by the impact, and construct the velocity and solid angle 
distributions of the stripped material. 

We find that the main sequence and subgiant companions lose 0.15 — 0.17 Mq as 
a result of the impact of the supernova shell, 15% of their mass. The red giant 
companions lose 0.53 — 0.54 Mq, 96% — 98% of their envelopes. The characteristic 
velocity of the stripped hydrogen is less than 10 3 km s _1 for all the scenarios: 420 — 590 
km s _1 for the red giant companions, 820 km s _1 for the main sequence companion, 
and 890 km s _1 for the subgiant companion. The stripped hydrogen and helium 
contaminate a wide solid angle behind the companion: 115° from the downstream axis 
for the red giant, 66° for the main sequence star, and 72° for the subgiant. With such 
low velocities, the bulk of the stripped hydrogen and helium is embedded within the 
low- velocity iron of the supernova ejecta. The hydrogen and helium may be visible in 
the late-time spectra as narrow emission lines. 

Although most of the stripped material is ejected at low velocities, all the numerical 
simulations yield a small high-velocity tail. The main sequence, subgiant, and the 
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red giant companions are just under the Delia Valle et al. (1996) upper limit from 
observations of SN 1990M taken near maximum light. 

The main sequence companion receives a kick of 86 km s™ 1 and the subgiant 
receives a kick of 49 km s _1 . In all cases, the kick to the remnant is smaller than the 
original orbital velocity. Because it is too small to intercept more than a negligible 
amount of momentum, the red giant core will not receive an appreciable kick. 

The impact of the supernova ejecta with the secondary star creates a hole in the 
supernova debris with an angular size of ~ 30° in the high-velocity ejecta and with 
an angular size of ~ 40° in the low-velocity ejecta. This corresponds to 7 — 12% of 
the ejecta's surface. Because we explore binary scenarios that are close enough, or 
almost close enough, to be in Roche lobe overflow, the degree of asymmetry is similar 
for all the models. The asymmetry in the supernova debris could have observational 
consequences beyond the change in morphology of the supernova remnant. The 
asymmetry in the supernova atmosphere will result in distorted P Cygni profiles that 
might indicate the presence of the companion star, but it will be difficult to use the 
degree of asymmetry alone to discriminate between a main sequence, subgiant, or red 
giant companion. 

The impact of the supernova shell will have consequences for the future evolution 
of the secondary star. After the impact, the main sequence star is puffed up, much like 
a pre-main sequence star. The luminosity will rise dramatically, to as high as ~ 5000 
Lq, as the extended envelope relaxes back into thermal equilibrium. The subgiant 
companion will follow a similar sequence of events. The star will not be contaminated 
by much supernova debris from the initial impact, but it may accrete low-velocity 
iron-group elements (or oxygen and silicon if the ejecta is radially mixed) at late times. 

A He pre-white dwarf will be left behind after a supernova explosion with a red 
giant companion. Almost all the envelope will be ejected by the impact, but a residual 
amount of material (~ 0.02 M©) will form an extended, hydrogen-rich envelope around 
the degenerate core. The star will evolve away from the red giant branch on a timescale 
of 10 5 — 10 6 years. On its post-RG track, it may appear as an underluminous O 
or B star before passing through an sdO or sdB phase on its way to a standard He 
white dwarf cooling track. This binary scenario could be a possible pathway for the 
formation of a subset of single, low-mass He white dwarfs. 

Subject headings: supernovae: general 
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1. Introduction 

Unlike a Type II explosion, a Type la explosion, the thermonuclear explosion of a CO white 
dwarf, is characterized by the absence of hydrogen in its spectrum. Although a massive star that 
is responsible for a Type II supernova explosion can be in a single or a binary system, the Type 
la progenitor scenarios that are based on the thermonuclear explosion of an accreting white dwarf 
require a companion as a mass donor to enable the white dwarf to reach its critical mass. However, 
despite searches for the remnant of the companion star, no direct evidence has yet been observed 
( Ruiz-Lap uente 1997] ). 

One method of discriminating between the many Type la progenitor scenarios is by searching 
for contaminating hydrogen and helium stripped from the companion star. However, this requires 
understanding the effect of the impact of the supernova blast on different companion stars to 
predict the amount of hydrogen stripped and its distribution in velocity and solid angle for the 
types of binary scenarios that have been proposed as progenitor models. To date, there have 
been only three numerical investigations of supernova impacts on companion stars. Fryxell & 
Arnett (1981) and Taam & Fryxell (1984) focused on impacts with low-mass main sequence stars. 
Livne et al. (1992) focused on red giant companions. In this paper, we improve upon the earlier 
theoretical work by exploring the nature of the interaction between the supernova ejecta and the 
companion star for specific Type la binary scenarios with higher resolution than was possible in 
the past. We select realistic ejecta profiles and realistic stellar models. It is no longer necessary to 
represent the supernova-secondary interaction by the collision of a simple shell with a polytropic 
companion. Although the earlier simulations studied the stripped mass and kick, they did not 
study or provide the distribution of stripped hydrogen in velocity and solid angle, which are 
essential to predict any spectral signatures of the impact. 

The goal of this project is to simulate with a 2-D hydrodynamics code the impact of a 
realistic supernova blast on a companion star in the context of current Type la progenitor 
scenarios in which the companion is hydrogen-rich. We select several of these scenarios from the 
literature to cover the possible range of secondary types which are typically low-mass (1.0 — 2.0 
M & ) main sequence, subgiant, or red giant stars. We use realistic stellar models and supernova 
debris profiles to represent each binary system. For each scenario we explore the hydrodynamics 
of the supernova-secondary interaction, calculate the amount of stellar material stripped from 
the secondary and the kick delivered by the impact, and construct the velocity and solid angle 
distributions of the stripped material. With the velocity and spatial distributions of the stripped 
material within the supernova debris it is possible to predict whether the contaminating hydrogen 
and helium are observable in the supernova spectrum as broad emission lines near maximum light 
or as narrow emission lines at late times. In addition, we can explore the change in structure of 
the secondary star caused by the impact to explore the possible final states (perhaps observable) 
of the secondary. 

We find that our main sequence and subgiant companions lose 0.15 — 0.17 Mq as a result 
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of the impact of the supernova ejecta. In contrast, the red giant companions lose almost their 
entire envelopes, 0.53 — 0.54 M@. The bulk of the stripped hydrogen, which has a characteristic 
velocity of less than or equal to 10 3 km s _1 , is embedded within the inner, low- velocity supernova 
debris. The complication that hydrogen and helium can be hidden deep and unseen within the 
ejecta, and obscured by competing iron and cobalt lines, can be addressed by radiative transfer 
calculations ( Pinto et al. 1999 ). In principle, radiative transfer calculations can reveal how much 
of the stripped material can be hidden and when it is likely to be observable. 



1.1. History of Theory 

The effect of a supernova explosion on a nearby companion star has been the subject of 
speculation for many years. One issue raised was the strength of the kick delivered to the 
companion. Colgate (1970) proposed that the shock-heating of the stellar envelope would cause 
material to evaporate preferentially back in the direction of the supernova. The secondary would 
receive a substantial kick from the ablation (evaporation from the surface of the star), as well 
as a smaller kick from the direct collision with the ejecta. Other groups QCheng 1974| ; |Wheeler| 



et al. 1975) made similar arguments supporting the suggestion that the kick due to ablation was 
larger than the kick from the collision. Numerical simulations by Fryxell & Arnett (1981) and 
Taam & Fryxell (1984) showed that the momentum transfer to the secondary was much less than 
previously thought and that ablation did not substantially increase the kick. The small binary 
separations we employ for the simulations in this paper, necessary for Roche lobe overflow, result 
in kicks larger than the ~ 10 km s -1 kicks reported by Fryxell & Arnett (1981). However, our 
typical 40 — 80 km s _1 kicks are still small in comparison with typical orbital velocities. 

The impact of the supernova debris on a nearby companion may be quite dramatic. The 
supernova ejecta may either directly strip material from the companion by direct transfer 
of momentum, or, through the conversion of the blast kinetic energy into internal heat, by 
evaporation. Because the incident energy (a fraction of 10 51 ergs) greatly exceeds the secondary's 
binding energy (10 48 ergs), Sofia (1967) suggested that a main sequence secondary in Roche lobe 
overflow would be completely destroyed. It was not until a more detailed study of the shock 
propagation into the stellar envelope was done that it was realized that, because the kinetic energy 
of the impact is deposited in the outer layers of the star, a main sequence star could easily survive 



such an impact ( Cheng 1974 ). 



Wheeler et al. (1975) analytically estimated the amount of mass stripped and the kick received 
by the companion as a result of the inelastic collision and the shock-heating. They calculated the 
outcome of the impact for an arbitrary secondary star using a geometrical parameter they devised, 
which depended on the incident momentum and the mass and compactness of the secondary. 
Their analysis implied that main sequence stars (represented by an n = 3 polytrope), being 
tightly bound, are only slightly stripped by an explosion in a close binary system. Red giants with 
loosely-bound envelopes are so catastrophically heated by the impact that the entire envelope is 
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ejected. 

The analytic prescription for the stripped mass and kick proposed by Wheeler et al. (1975) 
was put to the test in a series of numerical simulations of supernova impacts on both low-mass 
mam sequence companions ( |Fryxell fe Arnett 198 1| ; |Taam fe Fryxell 1984Q and low -mass red 



giants ( Livne et al. 1992] ) . Fryxell & Arnett (1981) demonstrated that the collision of a supernova 



shell with a 2.0 Mq, n=3 polytrope, with a binary separation 5.9 times the stellar radius, ejects 
0.013 — 0.052 Mq from the secondary, a mass loss roughly consistent with the analytic estimates of 
Wheeler et al. (1975). Although in their study the emphasis was on the momentum transfer rather 
than the stripped mass, Taam & Fryxell (1984) performed numerical simulations with n= 3/2 
and n= 3 polytropes that were designed to represent the range of compactness expected from a 
fully convective star to a radiative star. They found that the momentum was more efficiently 
transferred for the n=3/2 star, which is less centrally concentrated than the n= 3 star. The 
stripped n= 3/2 star presents a larger geometrical area to intercept the momentum, even when a 
higher fraction of mass is stripped, so that the momentum transfer is more efficient. More recently, 
Livne et al. (1992) simulated a Type la supernova explosion of a white dwarf in a close binary 
with a low-mass red giant that is almost close enough to be in Roche lobe overflow. The blast 
stripped the red giant of almost its entire envelope, imparting a velocity to the stripped material 
well below the velocity of the supernova ejecta. This not only confirmed the conclusion of Wheeler 
et al. (1975) that the impact would destroy the red giant, but also verified a prediction of Chugai 
(1986). Chugai (1986) focused on the supernova red giant interaction for Type la events in which 
the hydrogen donor is a symbiotic star - with a binary separation 3 — 10 times the red giant's 
radius. Chugai's analytic model predicted that 0.3-1.0 M© would be ejected from the companion 
with a characteristic velocity less than or equal to 10 3 km s _1 , much smaller than the 10 4 km 
s _1 characteristic velocity of the supernova ejecta, and that it would fill the inner 20% of the 
radius of the supernova ejecta. The confirmation by Livne et al. (1992) of the low velocity of the 
stripped hydrogen has important implications for the identification of Type la progenitors. Our 
more detailed calculations build on and extend this earlier, seminal work. 



1.2. Discriminating Between Type la Progenitor Scenarios by Searching for 

Hydrogen 

The exact amount of hydrogen stripped and its characteristic velocity is a key issue in light of 
the current interest in identifying the progenitors of Type la supernovae. Because it is difficult to 
reconcile the extremely hydrogen-poor spectra with the presence of a hydrogen-rich companion, the 
amount of stripped hydrogen can be used to discriminate between Type la progenitor scenarios. 
Type la explosions, which are characterized by a lack of hydrogen, could be the result of mass 
accretion onto a white dwarf in a binary system (single-degenerate scenario) flWhelan Sz lb en 1973| ; 
Nomoto 1982a) or the merging of two degenerate white dwarf stars (double-degenerate scenario) 



( Iben fc Tutukov 1984 ; Webbink 1984 ). Naturally, the double-degenerate scenario could account 
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for the lack of hydrogen in the supernova spectrum. However, no double-degenerate binaries 
have been found that are close enough to merge within a Hubble time and are massive enough 
to exceed the Chandrasekhar mass ([Marsh et al. 1995| ; |Saffer et ^.7998] ). The single-degenerate 



scenario is currently in favor ( [Livio 1999| ), If the companion is a main sequence, subgiant, or red 
giant star, the single-degenerate scenario necessarily implies that the binary system is rich in 
hydrogen. Hydrogen-rich material could be ejected from the secondary star as a consequence of 
the impact of the supernova shell. Hydrogen could also be present in the immediate environment 
as circumstellar material from stellar winds, mass lost from the primary in an earlier phase of mass 
transfer, or even as a layer of hydrogen on the white dwarf primary. If the hydrogen is swept up 
by the supernova ejecta, the origin of the material may determine its characteristic velocity, which 
in turn determines when it is most likely to be observed. Circumstellar material swept up with 
the supernova ejecta is expected to be observed as transient narrow H a emission or absorption 
lines ( [Wheeler 1992| ; [Filippenko 19971 ) near maximum light. If the hydrogen is stripped from the 



secondary and embedded within the inner iron layer of the supernova ejecta, it is more likely to 
be observed in emission as narrow H Q lines months after maximum light ( Chugai 1986j ), Searches 



for hydrogen must be targeted at either early-times near maximum light or at late-times when the 
photosphere has receded to reveal the iron layer (|Chugai 1986|) . 



Identification of hydrogen at early-times, either in emission or absorption, has been claimed 
in at least three Type la supernovae: 1981B, 1990M, and 1990N. Branch et al. (1983) presented 
high-quality optical spectra of SN1981B from maximum light to 116 days post-maximum. Their 
March 13 observation showed a small, narrow emission feature consistent with the rest wavelength 
of H Q . The feature was absent from a spectrum that was taken just 5 days later. Unable to 
firmly conclude that the feature was H Q , they suggested that high-resolution spectra be taken 
of supernovae near the same phase (~ 6 days past maximum light) to search for H a emission. 
Cumming et al. (1996) searched for narrow Hq, emission in high-resolution spectra of SNIa 1994D 
10 days before maximum and 6.5 days after maximum without detecting any hydrogen either 
in emission or absorption. Likewise, Ho & Filippenko (1995) in their high-resolution echelle 
observations of SNIa 1994D at 23 days past maximum did not detect any H Q in emission or 
absorption. 

From a series of observations of SNIa 1990N from 14 days before maximum light to 1 week 
past maximum light, Leibundgut et al. (1991) suggested that an unidentified absorption feature 
at 6300A could be hydrogen. They suggested that a thin layer of hydrogen on the surface of 
the progenitor star could be responsible for the feature, but because the velocity of the feature 
(1.2 x 10 4 km s _1 ) was much lower than that of the other lines, this was unlikely. Later, based on 
an analysis of optical and UV spectra of SN 1990N using synthetic spectra, Jeffery et al. (1992) 
concluded that the unidentified absorption feature at 6300A, and an additional unidentified feature 
at 6900A, were probably due to C II at 6580A and 7243A . 

Polcaro & Viotti (1991) identified a broad H a absorption feature in SNIa 1990M at four days 
after maximum light. This was later convincingly ruled out by Delia Valle et al. (1996), who 
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showed from a spectrum obtained a few days past maximum light that the absorption feature was 
due solely to observational bias. They found no evidence of H a in SN 1990M. Using an upper 
limit to the H a equivalent width and a measurement of the Si II 6350A equivalent width, they set 
an upper limit in the supernova atmosphere for H/Si of 2.0 x 1CP 6 , relative to solar. By assuming 



perfect mixing and using the abundances of the standard W7 model (Nomoto et al. 1984), they 
derived the upper limit to the mass fraction, X(H), of 0.007. However, they pointed out that if the 
hydrogen is buried deep within the supernova ejecta instead of being perfectly mixed, the integral 
mass fraction could be much higher. Likewise, if the hydrogen is restricted to the outer layers, the 
mass fraction could be much lower. Scaling X(H) by the ratio of the mass in the atmosphere to 
the total mass and assuming that 20 days past the explosion only 1/20 of the ejecta are revealed, 
they derived an X(H) of ~ 0.007 M atm /M tot ps 3 x 10~ 4 . 

A similar argument was made earlier by Applegate & Terman (1989) to place an upper limit 
on the hydrogen mass fraction in SN 1981B. Using the upper limit of [H/Si] < —4.0 in the LTE 
spectral fit to the maximum light observations of SN 1981B of Branch et al. (1982), Applegate 
& Terman (1989) derived an upper limit on the hydrogen mass fraction near maximum light in 
SN 1981B of 0.02. To find a lower limit to the hydrogen contamination due to the impact of 
the supernova ejecta on the secondary star, Applegate & Terman used the analytic prescription 
developed by Wheeler et al. (1975). They focused on a cataclysmic variable system with a 0.2 
Mq main sequence companion, represented by an n=3/2 polytrope, that was close enough to 
the primary to be in Roche lobe overflow. If the Type la progenitors are cataclysmic variables, 
Applegate Terman 1989) concluded that X(H) > 0.01. This is just below the Applegate &i 



Terman upper limit of X(H) < 0.02. In contrast, if the companion is a low-mass red giant 
with an envelope mass of 0.5 M© ; then X(H) > 0.20, which clearly would exceed the upper 
limit they deduced from observations of SN 1981B. Applegate & Terman (1989) conclude that 
binary scenarios with hydrogen-rich companions are not likely candidates for Type la progenitors. 
However, to directly compare their lower limit with their upper limit based on maximum light 
observations of SN 1981B, they implicitly assumed that the stripped hydrogen from either the 
low-mass main sequence star or the red giant companion was ejected at high velocities. 

It is currently believed, based on the work of Chugai (1986) and Livne et al. (1992), that 
any material stripped from a red giant companion will have a characteristic velocity <^ 10 3 km 
s , much lower than that of the supernova ejecta. As we show below, we agree with this general 
conclusion. If any hydrogen is embedded within the supernova ejecta, narrow hydrogen emission 
lines will appear only after the photosphere retreats into the iron layer. Thus, a critical test of 
the presence of a secondary star is the detection of low- velocity hydrogen lines in the late-time 
spectrum ( Ruiz-Lapuente et al. 199*3 ). However, such a detection is difficult because of the number 



of Fe lines that overwhelm the spectrum. To date, there has been only one claim of a detection of 
a low-velocity hydrogen line, ostensibly from a stripped companion. Ruiz-Lapuente et al. (1993) 
initially identified a weak H a line in a spectrum of SNIa 1991bg, a peculiar underluminous Type la 
supernova, 197 days after maximum light. However, this was not verified by Turatto et al. (1996). 
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Turatto et al. were not able to associate any emission lines with H a in their observations of SN 
1991bg, although they suggested a possible blend with [Co III] at 6578A or numerous Fe II and 
[Fe II] lines in that region of the spectrum. Garnavich & Challis (1997), in a reanalysis of a 
spectrum of SN 1991bg at 200 days after maximum light, confirmed the existence of the the 
narrow emission lines of Ruiz-Lapuente et al. (1993), but they could not unambiguously identify 
them. 

A key question for theoreticians is how much contaminating hydrogen and helium, whether of 
circumstellar or companion origin, can be hidden within the supernova ejecta. Very few theoretical 
upper limits have been set ( |Wheeler fc Harkness 19~9C ; Wheeler 1992j ). Wheeler (1992) claimed 



that, if the hydrogen is in LTE in the outer layers of the supernova ejecta, as much as 0.1 Mq could 
be present without contributing to the spectrum at maximum light. Branch et al. (1991) found a 
tighter upper limit of 1% by mass (~ 0.01 Mq) in a non-LTE calculation. By adding hydrogen to 
a synthetic spectrum and comparing with the maximum-light optical spectrum of SN 1981B, they 
concluded that 1% by mass is an upper limit on the hydrogen contamination for a layer of hydrogen 
uniformly mixed in the line- forming layer at maximum light. They found that, not only would 
a weak H Q line be visible, but, as an indirect consequence of hydrogen contamination, that the 
optical depth in the Lyman continuum would change the ultraviolet radiation field, the ionization 
structure in the line-forming layer, and, hence, the strengths of many of the non-hydrogen lines. 
Because Applegate & Terman (1989) estimated X(H) > 0.01, Branch et al. (1991) concluded that 
accreting white dwarfs with hydrogen-rich companions are unlikely to be progenitors of Type la 
supernovae. However, there are several caveats. First, the synthetic spectra employed for these 
investigations were done only near maximum light, when only high- velocity hydrogen would be 
visible. If most of the contaminating hydrogen is at low velocities, synthetic spectra need to be 
calculated at much later times when the photosphere has receded into the inner ejecta. Second, 
the hydrogen buried within the supernova ejecta may not be uniformly mixed in angle. Because 
non-uniform mixing may act to hide the embedded hydrogen, the angular dependence of the 
spectra needs to be considered. 

In this paper, we explore in detail the theoretical expectations for the distribution of hydrogen 
and helium in Type la debris, the hydrodynamic character of the impact, the kick to the secondary, 
and the nature of the post-impact structure. In §||, we review the numerical methods used. In 
§||, we describe the classes of supernova Type la progenitor scenarios and the criteria for selecting 
the candidate scenarios. In §||, we discuss details of the main sequence simulations, including 
the hydrodynamics of the impact, the stripped mass, the velocity and solid angle distributions 
of the stripped material, and the kick received by the companion. In §|5|, we present a sequence 
of simulations with the main sequence secondary to determine systematic trends in the quantity 
of stripped hydrogen and its velocity and solid angle distributions as the binary separation 
increases. In §||, we discuss details of the subgiant simulations, including a description of the 
impact, the stripped mass, the velocity and solid angle distributions of the stripped material, 
and the kick received by the companion. Section ^ contains similar details for the red giant 
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simulations. In §||, we explore the implications of the revealed velocity distributions of the 
stripped material and the possibility that contaminating hydrogen may be observable in a Type 
la supernova spectrum. In §|9|, we speculate on the effect of the impact on the future evolution of 
the secondary, and in §[M] we summarize our conclusions and suggest directions for future work. 
Postscript images, MPEG movies, and a selection of figures presented in this paper, are posted at 



http: / /www.astrophysics.arizona.edu and are available via FTP at www.astrophysics.arizona.edu, 



in directory pub/marietta. 



2. Numerical Methods and Techniques 

The hydrodynamics code we employ for these simulations is an extension of the code, 
Prometheus ( Fryxell et al. 198^ ), which is based on the Piecewise-Parabolic Method (PPM) 



flColella fc Woodward 1984j ; |Colella fe Glaz 1985|) . O ur version of PPM is a non-relativistic, 



explicit, automatically conservative, Eulerian scheme that achieves second-order spatial and 
temporal accuracy. Fluxes at interfaces are obtained by solving the Riemann shock-tube problem 
and shocks are resolved to one or two zones. 

A simple equation of state, including just radiation and ideal gas, is employed for all the 
calculations involving the main sequence companion and the envelope of the red giant; the 
degenerate core of the red giant is added only as a gravitational point mass. For the calculations 
involving subgiant companions, where degenerate electrons are required, we use a tabulated 
equation of state with arbitrary degeneracy and relativity. We tabulate the pressure, energy, 
and entropy of the ionization electrons, and any pair-produced electrons and positrons, by 
directly integrating the appropriate Fermi-Dirac integrals, following the prescription of Cox &; 
Giuli (1968). Using the ionization electron density (n e = pY e N av , where p is the density, Y e is 
the electron fraction, and N av is Avogadro's number) and the temperature (T) as independent 
variables, we use an iterative technique to solve for the degeneracy factor (r/ = fi/ksT, where 
[i is the chemical potential, and ks is Boltzmann's constant), which can be used to directly 
integrate the Fermi-Dirac integrals for pressure and energy. This method, although laborious, 
is thermodynamically consistent. The table, once generated, is accessed by a second-order 
interpolation routine. We estimate the pressures and energies to be accurate to better than 
~ 10" 4 . 

We incorporate into PPM an integral Poisson solver created by Miiller & Steinmetz (1995) 
to calculate the gravitational potential for an arbitrary 2-D mass distribution. Having found 
the gravitational potential, the gravitational force can be calculated using a finite difference 
approximation to the gradient. In the simulations involving the red giant companions, the 
gravitational point mass of the compact core is added directly to the force. In our implementation 
the solver expands the gravitational potential in Legendre polynomials. Up to twenty moments 
can be employed, including the dipole moment which we include to allow the secondary to 
move freely down the hydrodynamic grid in response to the impact of the blast. With the 
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addition of an interpolation routine, the hydrodynamical equations can be differenced in spherical 
(r,9), Cartesian (x,y), or cylindrical coordinates (p,z) while the gravitational potential can be 
constructed in spherical coordinates (r, 9). To improve accuracy, the gravitational potential can be 
centered on the center of mass at each cycle, in effect following the star as it moves down the grid. 
This decreases errors, notably in the stellar core, that over the long timescale of the simulation can 
otherwise slowly decrease the star's terminal velocity. Even with this improvement, the terminal 
velocity slowly decreases at an average rate of ~ 2.8% per stellar sound crossing time in a typical 



simulation with a low-mass main sequence companion. Thus, as will be seen in §4^3, we determine 
the kick after the terminal velocity is reached (4000 — 5000 seconds), but before the small errors 
have time to grow appreciably. 

For simplicity, we use all the mass on the hydrodynamic grid, stellar and supernova, to 
calculate the gravitational potential. Because the grid is so large, typically 6R (R = the stellar 
radius) in the p direction, and, since the explosion is ~ 1R away from the upper boundary, about 
42% of the supernova debris passes through the grid. Only a small fraction of the ejecta directly 
collides with the secondary. The characteristic velocity of the supernova ejecta is so fast that the 
error to the kick from including the supernova material on the grid is only ~ 4 km s _1 , which is 
an error of ~ 5%. 

For this project, a selection of low-mass stellar models were evolved by Chaboyer (1998). 
Following the lead of Sills et al. (1997), we reintegrate the main sequence and subgiant secondaries 
using a 1-D fourth-order Runge-Kutta program, the entropy and composition profiles, and the 
same equation of state used in our PPM code. We do not attempt to reproduce the distortions in 
the secondary due to the Roche lobe geometry. 

The weakly-bound envelope of the red giant requires special treatment. We reintegrate its 
entropy profile with a softened potential (cp(r) = —GM c /(r + r c )) to help stabilize it. Even with 
this adaptation, we find it necessary to switch to spherical coordinates for the red giant impacts. 
We surround each secondary star with a high-entropy, very low-density hydrostatic envelope 
which serves to fill the Eulerian grid with a background medium. In later sections we refer to this 
background material as the "circumstellar" medium to distinguish it from the stellar envelope. 
The new 1-D models with their hydrostatic envelopes are interpolated onto a 2-D cylindrical grid 
(main sequence or subgiant models) or a 2-D spherical grid (red giant envelopes). We use 2-D 
cylindrical coordinates whenever possible in order to efficiently follow the stripped material. To 
stabilize the envelope of the red giant, we must use 2-D spherical coordinates. 

We verify that each secondary star remains in hydrostatic equilibrium by running a simulation 
without the supernova explosion. The main sequence and subgiant simulations were run for three 
sound crossing times. Because the envelope of the red giant will be completely disrupted in less 
than one sound crossing time, the red giant simulation was run for only one sound crossing time. 
We estimate the fractional change in radius by the fractional change in gravitational energy for 
the main sequence and subgiant secondaries. We find that the radius changes by at most ~ 1%. 
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In the red giant case, the radius changes by at most ~ 8%. 

We perform all of the simulations with the main sequence and subgiant secondaries in 
2-D using cylindrical coordinates (p, z), with the z-axis defined to be the direction joining the 
primary and secondary of the original binary system. We estimate that the error in neglecting the 
orbital motion of the secondary, necessary for a 2-D calculation, is ~ 10~ 2 , the ratio of a typical 
orbital velocity of 10 2 km s _1 to a characteristic velocity of the supernova ejecta of 10 4 km s _1 . 
The thermonuclear explosion of the white dwarf occurs exterior to the grid and the supernova 
ejecta, specified by a realistic density and velocity profile, flows onto the grid via time-dependent 
boundary conditions in assumed spherical homologous flow. The other exterior boundaries are 
subject to the constraint that material is allowed only to leave the grid. The interior boundary 
(p = 0) is reflecting. 

We switch to spherical coordinates for the red giant companions because only in these 
coordinates is the loosely-bound envelope of the red giant stable. We position the red giant 
envelope at the origin and add a gravitational point mass to the force to represent the degenerate 
core. We assume throughout the simulation that the degenerate core does not move as a result of 
the impact. We justify this by noting that, because of the extremely small solid angle subtended 
by the degenerate core and its very high areal density, it can receive only a negligible (<C 1 
km s _1 ) kick. The supernova explosion, which in this case occurs on the grid, is interpolated 
onto the grid from a post-explosion ejecta profile. To avoid numerical problems associated with 
high-Mach flows, we assume an initial temperature in the supernova ejecta high enough that the 
internal energy is ~ 7.7% of the specific kinetic energy.^ To minimize numerical problems at the 
center of the supernova where the expansion leaves a low-density interior, we employ a minimum 
temperature which decreases linearly with time to mimic the cooling of the interior. Outside 
of the supernova ejecta at the interface between the supernova ejecta and the "circumstellar" 
medium, a minimum temperature of 100 K is enforced. These numerical problems appear only 
in the spherical calculations. The cylindrical calculations in which the supernova is added by 
time-dependent boundary conditions are always well-behaved. 

With the change to spherical coordinates, we alter the boundary conditions. The outer 
boundary is straightforward. As in cylindrical coordinates, the outer radial boundary is 
transmitting, but subject to the constraint that material is allowed only to leave the grid. The 
inner radial boundary is fixed at a small, but non-zero, radius. This helps to stabilize the 
interior of the envelope and, in addition, helps to avoid Courant problems. The inner boundary 
is non-transmitting (no mass flux allowed) with a zero-velocity condition. The inner boundary 
seems to affect the flow only after the envelope has been stripped. After the mass-stripping phase 
of the interaction, we regrid the calculation by removing the innermost zones, usually the first 25, 
in effect moving the inner boundary outward about a factor of 10 in radius, which increases the 



3 Our initial temperature is a factor of 1.0 — 4.0 higher than the temperature provided in the supernova ejecta 
model hedtbll (|Woosley fc Weaver 1994]). 
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timestep by the same factor, allowing us to continue the simulations efficiently. The mass enclosed 
in the innermost zones is added to the gravitational point mass. The region removed is typically 
only 1% in radius, and therefore, only 10~ 4 % in computational volume. The boundary along 
the axis of symmetry is reflecting. For high-velocity flow, the axis of symmetry manifests minor 
numerical artifacts which, however, do not affect the conclusions of the work. 

In addition to determining the stripped mass and the kick given to the companion, we find 
the distribution of the stripped material in both velocity and solid angle. The velocity distribution 
can in principle be used with a radiative transfer code to predict when, and if, the contaminating 
hydrogen is observable in the supernova Type la spectrum and to place useful constraints on 
the type of binary scenarios likely to be responsible for Type la supernovae. The solid angle 
distribution can be used to determine by how much the stellar material lags behind the supernova 
ejecta and the size of the solid angle contaminated by the stellar hydrogen and helium. It may 
also be relevant for estimates of the polarization of the emergent light. As stripped stellar material 
flows through the outer boundary, we record its mass, velocity, entropy, and composition. We 
advect the electron fraction Y e and ion fraction Yj along with ten general composition labels, five 
of which are reserved for the Type I supernova ejecta ("hydrogen", "helium", "oxygen-group", 
"silicon", and "iron-group"), four for the companion star ("hydrogen", "helium", "oxygen-group", 
and "silicon/iron-group"), and one for the "circumstellar" medium. The Yj and Y e are required 
by the equation of state, but we advect the composition labels as mass fractions so we can keep 
track of the origin of the stripped companion material and can distinguish it from the supernova 
and "circumstellar" material with which it is intermixed. To find the velocity and solid angle 
distributions, we include the distribution of the stripped companion material on the grid, as well 
as the stripped material that has left the grid. The stripped mass left on the grid is usually 
extremely small, given the length of a typical simulation. 

For the main sequence and subgiant simulations we generally employ a ~ 300(p) x 600(2:) 
cylindrical grid, and for the red giant a ~ 600(r) x 300(#) spherical grid, with nonuniform zoning 
that is finest in the region surrounding the secondary and coarsest near the edge of the grid. To 
cover 6 stellar radii around the secondary and still resolve the center of the star well enough 
to maintain hydrostatic equilibrium, nonuniform zoning is necessary. This nonuniform zoning, 
although efficient, creates difficulties in quantifying the effect of increasing or decreasing the 
resolution. 

To test the effect of the resolution we employ a cylindrical grid that has uniform zoning in 
a rectangular region enclosing the secondary. Outside this region, the zones are slowly increased 
in size to keep the calculation tractable. For the main sequence calculation, we performed two 
simulations: one with a spacing of 8.0 x 10 3 km requiring a 300(p) x 550(z) grid and one with a 
slightly smaller spacing of 6.0 x 10 3 km requiring a 320(p) x 845(z) grid. We can not vary the 
"resolution" by a large factor because greatly decreasing the resolution destabilizes the secondary 
while greatly increasing the resolution creates enormous Courant and data storage problems. 
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Nevertheless, we find that changing the resolution in the region around the secondary can 
affect the results in several systematic ways. First, the amount of mass stripped increases, but 
only very slightly, as the resolution increases. In the high-resolution simulation the secondary 
lost almost ~ 0.1% more mass than in the low-resolution simulation. Second, as the resolution 
increases, the kick to the remainder of the secondary increases. The velocity of the center of 
mass follows the same overall profile, the only difference being a slight increase of 1.1% in the 
terminal velocity. An increase in momentum transfer with higher resolution was noted by Taam 
& Fryxell (1984). They attributed this trend to a more accurate description of the momentum 
transfer in the outer layers of the star, where the shock energy is deposited, because of the steep 
density gradient there. Third, as the resolution increases, the velocity of the stripped material 
systematically shifts to slightly lower velocities. The velocity at the half-mass point shifts from 
841 km s _1 to 829 km s _1 . If with higher resolution the transfer of momentum to the remainder of 
the secondary is more efficient, then less momentum is available for the stripped material. Apart 
from the subtle systematic shift as the resolution changes, the features in the moderate-velocity 
region (1000 — 3000 km s _1 ), which we associate in §4.1 with the composition transitions in the 
supernova ejecta, vary in position and amplitude. However, the velocity distributions of the 
low- and high-resolution calculations show the same overall profile, especially at the low- and 
high- velocity tails. Above 3 x 10 3 km s _1 , the low-resolution distribution has 12% more mass than 
the high-resolution simulation. To reprise, as the resolution increases, the amount of stripped 
mass and the kick imparted to the secondary increase slightly, and the half-mass velocity of the 
stripped material decreases slightly. 



3. Standard Type la Supernova Binary Scenarios 

A Type la supernova could be the product of a merger of CO white dwarf with a He or CO 
white dwarf (double-degenerate scenario), or the explosion of a white dwarf which has reached 
its Chandrasekhar mass by accreting hydrogen or helium from a nondegenerate companion 
(single-degenerate scenario) ( pben &: Tutukov 1984 ). These two competing binary scenarios 



originate from very different evolutionary paths and naturally have different predictions for the 
Type la supernova rates. Binary population synthesis, a statistical tool for exploring the possible 
phases of binary evolution, can be used to estimate rates for progenitor scenarios to discriminate 
between progenitor models (cf. [ben fc Tutukov 1984 [Tutukov et al. 1992j ; |Branch et al. 19951 ). 



However, the Type la supernova rates from population synthesis calculations suffer from serious 
uncertainties. Common envelope evolution is not well-understood. It is parameterized by acE, 
the efficiency with which the orbital energy is used to eject the common envelope, whose value is 
uncertain. The realization frequencies are dependent on the stellar population and age. Another 
complication arises from the different models for the thermonuclear explosion of the white dwarf. 
In a Chandrasekhar explosion (the standard model), the CO white dwarf accretes material until it 



approaches the Chandrasekhar mass and carbon ignites in the core (Wbosley & Weaver 1986). In 



sub-Chandrasekhar explosions, the ignition occurs in a helium layer around the CO white dwarf 
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before the Chandrasekhar mass is reached ( [Nomoto et al. 1982b ). In general, binary population 



synthesis calculations favor the double-degenerate scenario over the single-degenerate scenario 
( Branch et al. 1995| ). Recent calculations by Branch et al. (1995) for a young population (10 8 



yr) find double-degenerate rates of ~ 1CP 3 yr -1 . This is closer to the current galactic Type 
la supernova rate of 4.0 x 10 -3 yr -1 than the rates of 10~ 4 — 1CP 6 yr -1 for each of the many 
possible single-degenerate scenarios. But, for an older population (10 10 yr) the double-degenerate 
rate drops to ~ 10 -4 yr" 1 which can be matched by the realization frequency for symbiotic 
systems, a single-degenerate scenario in which the white dwarf accretes hydrogen from the wind 
of its red giant companion. Likewise, Ruiz-Lapuente (1996) finds that for an Sb galaxy the 
double-degenerate rate is ~ 10~ 4 yr -1 , but that several single-degenerate scenarios have rates 
similar to this. Because of the complexity of the population synthesis calculations, it is difficult to 
exclude any Type la progenitor scenario outright just on the basis of realization frequencies. The 
presence of hydrogen in Type la supernova spectra is a stronger indication of a single-degenerate 
progenitor system because any Type la explosion in a double-degenerate binary system is expected 
to be hydrogen- free. In this paper, we examine only the single-degenerate scenarios for Type la 
supernova explosions. 

Single-degenerate Type la models can be coarsely subdivided (see Table [j]) into Hydrogen 
Cataclysmic Variables (H CV), Hydrogen Cataclysmic-Like Variables (H CVL), Symbiotic Stars 
(SS), Hydrogen Algols (H Algols), Helium Cataclysmic Variables (He CV), and Helium Algols 
(He Algols), based on the evolutionary stage of the secondary, method of mass transfer, and 
composition of the secondary's envelope ( Branch et al. 1995| ; [Ruiz-Lapuente 1996| ). The secondary 



can be a low-mass main sequence, subgiant, or red giant star that slowly loses its envelope to the 
primary by Roche lobe overflow for close binaries or by stellar winds for wider binaries. The mass 
transferred can be hydrogen or helium in cases where the hydrogen layer has already been lost by 
earlier phases of mass transfer. 

In the H CV scenario, a CO white dwarf, left in a close binary orbit by an earlier episode 
of common envelope evolution in its AGB phase, accretes hydrogen by Roche lobe overflow from 
a lower-mass main sequence secondary. The mass transfer is maintained by magnetic braking 
and has a characteristic rate of M > 1.0 x 10~ 8 Mq yr -1 ( Verbunt &: Zwaan 198l| ). Not all 



H CV systems end as Type la supernova progenitors. Only the largest white dwarfs with the 
highest mass transfer rates, high enough to suppress hydrogen flashes, can accrete enough mass 
to reach the Chandrasekhar mass. Livio & Truran (1992) associated these binaries with recurrent 
nova systems. When the Chandrasekhar mass is reached, the white dwarf explodes as a Type la 
supernova, but only if it can avoid accretion-induced collapse ( Nomoto fc Kondo 199C)| ). Delia 



Valle & Livio (1994) suggested these binaries could be responsible for Type la events in late- type 
galaxies while another progenitor from an older stellar population could be responsible for Type 
la events in early-type galaxies. 

The H CVL systems have an evolutionary path like that of the H CV systems. After an 
episode of common envelope evolution, a white dwarf is left in a close binary orbit with a slightly 
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higher-mass main sequence or subgiant secondary. Because of the mass ratio, the mass transfer 
is unstable and proceeds on a thermal timescale. This higher mass transfer rate (M ~ 10~ 6 Mq 
yr -1 ) is responsible for the steady burning of hydrogen to helium on the surface of the white dwarf. 



Eventually, the Chandrasekhar mass is reached and the star explodes (Rappaport et al. 1994 
Hachisu et al. 1996; see also [Hachisu et al. 1999b ). 



We use the H Algol and Symbiotic classification, following the path of Branch et al. (1995), 
as a general term to describe binaries in which the donor is a red giant that slowly loses mass to 
the white dwarf. In these scenarios, after a common envelope phase, a CO white dwarf is left in 
a binary orbit with a main sequence star. In the H Algol scenario, the main sequence secondary 



evolves all the way to the red giant phase before filling its Roche lobe for the first time ( Whelan 



& Iben 1973; Iben & Tutukov 1984). Mass transfer to the white dwarf is by Roche lobe overflow, 
like the other Type la progenitor scenarios. In the symbiotic scenario, the red giant is too distant 
to fill its Roche lobe in the course of its evolution. Mass loss will occur by stellar winds. Many 
evolutionary paths can lead to the symbiotic scenario, the simplest being a wide binary whose stars 
evolve independently. Other possible paths include prior phases of common envelope evolution or 
conservative mass transfer during the primary's evolution, which leaves the binary separation too 
wide for the secondary to ever fill its Roche lobe ( Yungelson et al. 1995| ). Until recently, symbiotics 



were considered as viable Type la candidates only for sub-Chandrasekhar explosions ( Kenyon 



et al. 1993| ; Yungelson et al. 1995| ). However, recently Hachisu et al. (1999a) have proposed that 



symbiotic systems with Chandrasekhar explosions can account for the Type la rate. 

The He CV and He Algol classes are analogous to the H CV and H Algol classes, the primary 
difference being that the mass transferred is helium instead of hydrogen. This complicated 
evolutionary path begins with the creation of a white dwarf in a close binary orbit, following a 
phase of common envelope evolution with a main sequence star, as in the previous scenarios. 
The secondary evolves though the red giant phase. However, unlike the earlier scenarios, the 
secondary catastrophically loses its hydrogen envelope in another common envelope phase. The 
secondary, now a helium star, continues to evolve. For the He CV scenarios, mass transfer begins 
when the helium star goes into Roche lobe overflow. Stable mass transfer (q < 1) is maintained 
by gravitational wave radiation and proceeds on a nuclear timescale. For these binaries, helium 
star masses are typically in the range of 0.3 — 1.0 Mq and white dwarf masses are in the range of 
0.6 — 1.0 Mq. These are Type la supernova Chandrasekhar and sub-Chandrasekhar explosions 



(Iben & Tutukov 1991). For the He Algol scenarios, mass transfer begins during the secondary's 
expansion to a helium giant. Mass transfer is quasi-stationary (M~ 10~ 6 Mq yr _1 ) with secondary 
masses in the range of 0.75 — 2.3 Mq and with white dwarf masses in the range of 0.9 — 1.4 Mq. 



The He Algol systems are candidates for Type la Chandrasekhar explosions (Iben & Tutukov 



1994). Branch et al. (1995) finds that the realization frequencies for these scenarios are 10" 4 -10 



yr 1 for young populations and are negligible for old populations. 

In this paper, we focus on exploring the supernova-secondary interaction for the four 
hydrogen-rich scenarios in Table To simulate each impact, we select a representative binary 
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system from the literature for each of the four scenarios and estimate the binary separation and 
the companion's mass at the time of explosion. We use a stellar model in the simulation that is 
reasonably close in mass and evolutionary stage. We use either SNIa W7 ( [Nomoto et al. 1984| ) to 
represent a Chandrasekhar explosion or SNIa Hedt ( Woosley fc Weaver 1994 hedtbll) to represent 



a sub-Chandrasekhar explosion, depending on which is more likely for each binary scenario. The 
secondary models are listed in Table ^ with the masses and radii provided. Information on the 
supernova ejecta models is given in Table |. 

To represent the H CV class, we construct a prototype binary, based on the Type la 
supernova progenitor candidates of Livio & Truran (1992), consisting of a 1.0 Mq main sequence 
secondary that lost 0.1 Mq to a massive 1.3 Mq white dwarf primary, allowing it to reach the 
Chandrasekhar mass. Assuming that the binary is in Roche lobe overflow at the time of explosion 
and that the secondary just fills its Roche lobe, from Eggleton's relation ( jEggleton 1983 ) for 



q = 1.0/1.4 we estimate a binary separation of ~ 3R. For this simulation, we represent the 
secondary by a 1.0 Mq solar model ( ptiaboyer 1998| ), the exploded white dwarf by the SNIa W7 



(Nomoto et al. 1984) ejecta profile, and employ a binary separation of 3R. 



We use the work of Li & van den Heuvel (1997) on supersoft X-ray sources as Type la 
supernova progenitor candidates to construct sample scenarios for both the H CVL and H Algol 
classes, one for each candidate proposed by Li & van den Heuvel (1997). Although H CVL 



binary systems have been included in the past as part of population synthesis studies ( Rappaport 



et al. 1994| ), they are now the subject of renewed interest due to the suggestion by Hachisu 
et al. (1996) that a white dwarf fed by an optically thick wind can experience a much higher mass 
transfer rate than previously thought. Hachisu et al. (1996) have suggested that these binary 
systems are the progenitors of some supersoft X-ray sources, as well as potential Type la supernova 
progenitors, although the realization frequency is in dispute (Yungelson & Livio 1998). 



As prototypes for our H CVL and H Algol classes, we select the two Li & van den Heuvel 
(1997) scenarios: the first scenario with a 2.5 Mq subgiant companion and the second with a 1.0 
Mq red giant companion, prior to mass transfer. We represent the subgiant companion just prior 
to mass transfer with a 2.1 Mq subgiant star having a radius of 2.3 Rq ( |Chaboyer 1998| ). To 
simulate the effect of binary mass transfer, which was not included when the star was evolved, we 
adjust the entropy profile to decrease the mass and radius of our secondary to 1.13 M Q and 1.7 
Rq , respectively, which reproduces the mass and radius of the secondary at the time of explosion 
as estimated from Li & van den Heuvel (1997). We do not attempt to reproduce the effect of 
the mass transfer on the red giant secondary. We represent the red giant companion with a 0.98 



Mq red giant with a radius of 170 Rq ( Chaboyer 1998 ). Because we can not include binary 



mass transfer in the stellar evolution code, our stellar models can not be exact recreations of the 
secondaries used in the binary evolution calculations of Li & van den Heuvel (1997). However, as 
we shall see, the stripped fraction and the velocity distribution of the stripped material are fairly 
robust and predictable. 
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To represent the Symbiotic (SS) class, we construct a prototype system based on the work 
of Yungelson et al. (1995), who concluded that symbiotic systems could be the progenitors of up 
to 1/3 of Type la events in young and intermediate age (t ^ 6 X 10 9 yr) stellar populations. We 
select the same 0.98 M Q red giant secondary as in the H Algol case, but to ensure that mass 
transfer takes place by stellar winds we increase the binary separation so that the red giant 
could not have filled its Roche lobe. Because these systems are good candidates for Type la 
sub-Chandrasekhar explosions, we select the ejecta model SNIa Hedt (hedtbll, |Woosley fc Weaver 



1994 ) to represent the exploded white dwarf. 



Table || lists the key simulations performed for this paper, with the stellar model, supernova 
ejecta model, and binary separation that were used in each simulation. The four simulations, 
representing the four hydrogen-rich classes of single-degenerate progenitor models, are labeled 
after their respective class as HCV, HCVL, HALGOLa, and SYMB. Additional calculations were 
performed (HCVa, HCVb, HCVc, and HCVd) to gauge the trends with separation in the mass 
stripped, the kick, and the velocity distribution of the stripped hydrogen. We include a simulation 
using a 2.1 Mq subgiant (HCVLa) secondary, which we compare to the HCVL simulation, whose 
secondary suffered substantial mass loss prior to the supernova explosion. 

Our choice of reference binary scenarios is motivated in part by the likelihood that a 
particular progenitor can be responsible for a significant fraction of Type la events, and in part 
by the desire to cover the range of parameter space likely to be encountered by any Type la 
explosion in a single-degenerate binary system. We expect that the four binary scenarios HCV, 
HCVL, HALGOLa, and SYMB will cover the basic range of secondary stars, supernova explosion 
models, and binary separations expected for Type la supernova progenitors with a minimum 
of overlap. For example, the HCV, HCVL, and the HALGOLa simulations represent low-mass 
main sequence, subgiant, and red giant companions in Roche lobe overflow, in order to cover the 
range of compactness in the the evolution of a low- mass star. The mass profile of each secondary 
model is shown in Figure |], and the binding energy profile is shown in Figure ^. The kinetic 
energy of the supernova blast is roughly 11, 10, 2300, and 1300 times the binding energy of the 
secondaries for the HCV, HCVL, HALGOLa, and SYMB simulations, the difference being the 
large difference in binding energy between the tightly-bound main sequence and subgiants and the 
loosely-bound envelope of the red giant. To explore the role of the supernova debris profiles, we 
employ two Type la models (see Table |||): W7, a standard Chandrasekhar explosion model, and 
Hedt, a sub-Chandrasekhar model with 78% of the momentum of W7. The HALGOLa and SYMB 
simulations represent the red giant companion with the same secondary model, but the HALGOLa 
simulation represents the supernova ejecta with W7. The SYMB simulation represents the ejecta 
with Hedt at a slightly greater binary separation. Because the mass transfer is wind-driven, the 
solid angle subtended by the red giant in the SYMB simulation is therefore slightly smaller. 
The solid angle subtended determines the fraction of the supernova momentum incident on the 
secondary. For any low-mass companion in Roche lobe overflow the fraction of the sky subtended 
is typically (l/4)(R/a) 2 ~ 1/36 because a~ 3R, independent of the companion's stellar type. 
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4. The Main Sequence Companion 

In this section, we describe the impact of the blast wave on a main sequence star, our HCV 
scenario. As in the subgiant and red giant cases, we neglect any changes to the structure of the 
star from the Roche potential. Likewise, we do not add an accretion disk, magnetic field, or any 
other intrinsic features of cataclysmic variables which would introduce complications. For this 
simulation, we use a 300(/>) x 575(z) cylindrical grid with the 1.0 M main sequence secondary, 
centered on the origin. Approximately 161 (p) x 322(z) zones are allocated to the secondary. 
The grid extends 6 stellar radii in the p direction, 12 stellar radii in the downstream (negative 
z) direction and 2 stellar radii in the upstream (positive z) direction. The white dwarf primary 
explodes just off the grid in the upstream direction, and the supernova debris flows through the 
upper boundary onto the grid by time-dependent boundary conditions. The density and velocity, 
which characterize the debris structure, energy, and momentum, are scaled with time: the Eulerian 
velocity as t~ l and the Lagrangian density as t~ 3 . 

The HCV was run for 2.0 x 10 4 seconds, slightly more than 7 of the companion's sound 
crossing times and much longer than the 1.4 x 10 3 seconds that it takes for the trailing edge of 
the supernova to sweep past the back of the companion. The length of the calculation is chosen 
to be sufficiently long, and the size of the grid sufficiently large, that we can follow the progress 
of the shock-heated material ejected from the companion in both velocity and solid angle and can 
construct a distribution for each. In addition, the long simulation time allows the companion 
star to begin to reestablish hydrostatic equilibrium. Although we can not follow the companion's 
Kelvin-Helmholtz evolution, its post-impact structure may provide us with some insight into the 
possible consequences of the impact for the companion's future evolution. 



For clarity, this section is subdivided into distinct topics. In §4.1, we describe the debris 



structure of the impact, and in §4.2 we consider the consequences of the impact for hydrogen 



burning. In § |4.3| , we present and describe the stripped mass for this scenario. In § |4.4j , we 
qualitatively discuss the possibility that the main sequence companion may accrete supernova 
material and contaminate its envelope with heavy metals. In § |4.5| , we discuss the velocity and solid 



angle distributions of the stripped companion material and possible observational implications. In 



4.6, we discuss the kick received by what remains of the companion. 



4.1. Hydrodynamics of the Impact on the Main Sequence Companion 

The impact of a blast wave on a main sequence star begins a sequence of events common to 
all supernova-companion interactions. The initial impact of the supernova shell on the surface of 
the star drives a shock into the stellar envelope. A reverse shock propagates back into the ejecta. 
A contact discontinuity marks the interface between the supernova ejecta and the shocked stellar 
material. The shock propagating through the envelope decelerates as it runs up the steep density 
gradient of the star. After the shock passes through the stellar core, it accelerates down the 
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decreasing gradient. The shock front, highly curved because of the density gradients, converges 
in the back of the star. The reverse shock develops into a bow shock around the companion star 
that smoothly deflects the rest of the incoming supernova debris around the companion. After 
most of the supernova debris has passed by, the outer layers of the stellar envelope, which have 
been shock-heated to such an extent that the new speed of sound exceeds the companion's escape 
velocity, are ejected, often corrugated by Kelvin-Helmholtz and Rayleigh- Taylor instabilities. 
The stripped material flows slowly away from the star, embedded within the inner layer of the 
supernova ejecta. The stellar core expands and cools in response to the initial compression by the 
shock. The star now has an extended, asymmetrical envelope and begins to pulsate to bring the 
envelope back into hydrostatic equilibrium. Figures ^ - [4] illustrate the sequence of the interaction 
in a series of 2-D images from the initial impact to the reestablishment of hydrostatic equilibrium. 
Figures || illustrate the mass stripping of the companion in a series of 2-D cartoon images, with 
each color indicating the dominant element in that region. 

Figure |3|a shows the companion centered at the origin of the cylindrical grid at the beginning 
of the simulation. Figure ^b shows the initial impact of the blast at ~ 100 seconds after the 
explosion. The supernova ejecta can be seen flowing onto the grid, sweeping up "circumstellar" 
material as the leading edge begins to flow around the secondary star. There is a region of high 
pressure at the site of the impact on the companion's surface where a shock is being driven into 
the star, and a bow shock is beginning to develop that will divert the flow of the supernova ejecta 
around the companion. The back of the companion is still undisturbed. 

The density variations in the ejecta that are visible as spherical rings in Figure ||b can be 
directly identified with the density variations in Figure ^ by their velocity. For example, the first 
strong ring striking the surface of the companion star in Figure ||b has a velocity of 1.59 x 10 4 
km s _1 and can be identified with the first sharp density spike in Figure || The second ring has 
a velocity of 1.36 x 10 4 km s _1 and can be identified with the second density spike in Figure 
^. These features are associated with changes in the composition of the ejecta: the transitions 
between the "oxygen", "silicon", and "iron" layers, (see Figure ||a). 

Figure ||c shows the impact 179 seconds after the explosion. The tail shock and bow shock 
dominate this image. By this time, the collision of the supernova ejecta with the "circumstellar" 
medium has created a curved shock front which has converged along the downstream axis after 
flowing around the companion. The shock front reflects at the axis and begins to move outward. 
The tail shock is visible in Figure ||c as the flared vertical discontinuity along the downstream 
axis. The upstream axis is dominated by the bow shock which wraps smoothly around the top of 
the companion, but becomes distorted farther down as a consequence of the impact of the density 
variations of the ejecta. The material between the bow shock and the shocked surface of the star 
is primarily supernova debris mixed with "circumstellar" medium. The pressure in this region now 
rivals the pressure in the interior of the star (see Figure [7]). 

At 329 seconds (Figure |3|d), the shock front passing through the companion has a strong 
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curvature caused by the density gradient of the star. The center of the shock decelerates as it 
encounters the rising density near the stellar core, while the wing of the shock accelerates through 
the back half of the companion where the density is decreasing. The shock-heated envelope of 
the companion starts to suffer mass loss; escaping material is now visible in the flow that passes 
around the star behind the bow shock. Kelvin-Helmholtz and Rayleigh- Taylor instabilities in the 
shocked envelope are just barely visible in this image. 

The ejected companion material is confined to a torus around the back of the companion - 
with the inside of the torus filled with "circumstellar" material that has been compressed behind 
the companion. There, the flow is almost directly parallel to the downstream axis (see Figure |5|c). 
This material is the first to be ejected from the star and is mixed with the oxygen layer of the 
supernova debris. This will form the high-velocity tail of the velocity distribution of the stripped 
companion material. In Figure |B|c, the iron shell has expanded to enclose the star, the silicon shell 
is just below the back of the star, and the oxygen shell is everywhere exterior to the silicon shell. 
More than 90% of the supernova's mass and momentum that will collide with the companion is 
now represented on the grid. 

Figures ||a and ||b show the interaction at 529 and 755 seconds after the explosion (see also 
Figures [5]d and ||e). By 529 seconds, the collision is effectively finished. Only a negligible mass and 
momentum will collide with the companion after this time. The shock within the companion has 
just passed through the stellar core and will begin to accelerate down the density gradient of the 
back of the companion. Kelvin-Helmholtz and Rayleigh- Taylor instabilities, the fingers of swirling 
material in the shocked envelope, are very visible by this time. As yet, only 0.07 M & has been 
ejected. By 755 seconds, the shock within the companion has converged along the downstream 
axis. In subsequent frames, a reflection shock can be seen progressing back through the remnant 
of the companion. 

Figure ^|c shows the end of the mass stripping phase. By 2027 seconds, the star has expanded 
and pushed the bow shock back towards the site of the explosion. About 93% of the companion's 
mass that will be ejected has been ejected by this time. The velocity field around the star is now 
more spherical. This can be inferred from Figure ||f by the change in the distribution of the stellar 
hydrogen, which is no longer confined along the downstream axis. About 1000 seconds later, the 
first in a series of shocks is driven into the envelope marking the beginning of a global stellar 
pulsation phase. 

Figure ||d shows the remnant after three oscillations, 2 x 10 4 seconds after the supernova 
explosion. The star is almost spherical again, but it is surrounded by a complicated system of 
shocks. The shock structures can be qualitatively explained in terms of the overtaking of each 
shock front by the subsequent shock front. Add to this mixture of shocks the reflected shock from 
the original shock through the companion, and the structure becomes difficult to disentangle. 

Exactly what happens to the structure of the companion as the shock is driven through the 
star? To reliably follow the companion from the impact to the reestablishment of hydrostatic 
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equilibrium, we ran a slightly different simulation, with the same binary scenario, but with a 
uniform grid and the gravity solver centered on the center of mass at each timestep, which is more 
suitable for long simulation times. Figures |7| - |9] show the evolution of the secondary's pressure 
profile from the beginning to the end of that simulation. Figure ^ shows the shock climbing the 
density gradient towards the stellar core. At ~ 400 seconds, the shock passes through the stellar 
core and begins to descend down the density gradient. Also, as a result of the impact, a reverse 
shock is formed at the interface between the star and supernova that propagates back into the 
supernova material. Eventually, the reverse shock becomes the bow shock. 

As Taam &; Fryxell (1984) discussed in their work on impacts on low-mass main sequence 
secondaries, the ram pressure of the ejecta is roughly equal to the central pressure of the main 
sequence secondary. The typical ram pressure of the ejecta is P = pv 2 M e jV 2 /(4-7ra 3 /3), (N.B., a 
~ 3R when it hits a low-mass secondary in Roche lobe overflow). For our main sequence simulation 
this is ~ 10 17 dynes cm -2 , very close to the central pressure of 2.37 x 10 17 dynes cm -2 . The 
ram pressure of the ejecta is close to the central pressure of the secondary only in the simulation 
with the main sequence companion. The more evolved subgiant and red giant have larger radii, 
which increases the binary separation and, thereby, decreases the incoming density and the ram 
pressure at first impact. For our subgiant and red giant secondaries, the expected ram pressures 
are roughly 10 16 and 10 10 dynes cm -2 , respectively, which is too low to significantly affect the 
subgiant's core by 2 orders of magnitude and the red giant's core by 7 orders of magnitude. In our 
HCV calculation, we find that the central pressure changes by a factor of only 2.5 as the shock 
passes through the core, which makes it a weak shock. The central density rises by only a factor 
of ~ 1.7, the temperature by ~ 1.5 and the entropy by 2.5 k^ 1 baryon -1 . 

After the initial pressure rise, the star expands and the central pressure begins to drop. The 
shock accelerates down the density gradient in the back of the star until at ~ 700 seconds it 
passes through the limb of the star and converges on the downstream axis. Figure || illustrates 
the dramatic pressure rise at the back of the star as the curved shock fronts converge. The high 
back-pressure exerts a force on the secondary which decreases the kick received from the impact. 
As part of the shock-heated envelope escapes from the star, the back-pressure drops, until after 
approximately 2000 seconds, the mass stripping phase is complete and the remainder of the 
secondary reaches its terminal velocity. At this stage, the star is slowly expanding. However, the 
passage of the shock and the ejection of part of the envelope begin a phase of radial pulsations 
during which shocks are sent out through the extended envelope. It is during this phase that the 
post-impact secondary recovers hydrostatic equilibrium. 

Figure |^ shows the pressure profile from the end of the mass stripping phase to the 
reestablishment of hydrostatic equilibrium. As slow moving material begins to settle back onto 
the star, a bounce occurs (~ 4000 seconds) and a shock is driven into the outer envelope. Every 
3000 — 4000 seconds, roughly a sound crossing time, another pulsation with a smaller amplitude 
is seen moving through the extended envelope. Despite the dramatic appearance of the shocks in 
the extended envelope, little or no mass ejection takes place as a result. As shown in Figure (TTJ, 
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the stellar core also pulsates, with a period of ~ 2600 seconds. In the last pressure profile shown 
in Figure |9[ the central pressure has stabilized. 

The changes in the central pressure, density, temperature, and entropy are listed in Table ||. 
All quantities but the entropy, which increases by shock-heating, have dropped significantly due 
to the expansion. Even after the interior of the star stops expanding and the infalling material 
compresses the core, the central pressure and density are not restored to their original values. 
Although such pulsations were not noted in earlier numerical simulations, probably due to shorter 
simulation times, Taam & Fryxell (1984) did find that some of their low-mass main sequence 
companions suffered an abrupt rise in central density followed by a severe drop as the star began 
to expand. The density and entropy profiles of the pre- and post-impact companion are shown in 
Figures |ll| and [l^. The radius of the remnant star, which we define as the radius that encloses 
about 90% of the remnant's mass, is approximately 7.4 x 10 5 km, a modest increase of less than 
10%. The remaining 10% of the mass forms a low-density extended envelope out to 6.7 x 10 6 km 
(~ 10 Rq) with a composition similar to that of the original stellar envelope. 



4.2. Quenching Hydrogen Burning in the Stellar Core 

Do the abrupt rise and fall of the central temperature and density have any unique 
consequences for the energy generation rate and the luminosity of the companion during the 
impact? We expect that the energy generation rate, with its extreme temperature sensitivity, 
should rise and fall dramatically as the stellar core heats and cools. Figure [l(] illustrates the rise 
and fall of the central temperature and density. The oscillations that the secondary experiences 
as hydrostatic equilibrium is reestablished are visible in both profiles. For a quick estimate of the 
energy generation rate, we calculate the hydrogen burning rate via the three pp chains and the 
CNO cycle flClayton 1985 ; Kippenhahn fc Weigert 1991 ) along a 1-D cut through the densest (and 



hottest) part of the companion. As can be seen in Figure H the energy generation rate shoots 
up by a factor of 1000 as a result of shock-heating, then drops precipitously as the star expands. 
Because the compression phase is so brief, despite the extremely high energy generation rate, only 
a small amount of additional energy is deposited in the center of the star. For the main sequence 
star an additional ~ 10 38 ergs has been deposited, which is extremely small in comparison with 
the binding energy of the star (~ 10 48 ergs). However, because the timescale of the compression 
phase is close to the sound crossing time of the stellar core, the star may be able to expand and 
cool. Hence, 10 38 ergs is an upper limit. By the end of the simulation, the energy generation rate 
is only 0.07 L . The hydrogen burning has been quenched by the impact. The implications of this 
are discussed in §|9[ 
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4.3. The Stripped Mass of the Main Sequence Companion 

The amount of mass stripped from the secondaries is shown in Table |6[ The third column 
lists the stripped mass for each scenario based on the numerical calculation and the fourth column 
lists the stripped mass based on the analytic estimation of Wheeler et al. (1975). To determine 
the amount of material ejected by the impact and the passage of the shock through the star, we 
compare the velocity of the companion's material within each zone with its local escape velocity 
y / 2^>, (ft being the gravitational potential. If its velocity exceeds the escape velocity, we consider 
the mass unbound. If not, we consider the mass bound to the star, only displaced, although, 
because of the Eulerian nature of the code, we are not able to follow the calculation until the 
displaced stellar material falls back onto the companion. When any stellar material reaches the 
edge of the grid, we assume that if it is unbound at that time it will remain unbound. By the end 
of the simulation, 0.15 Mq of the main sequence stars's envelope is unbound, and all but 2.4% of 
this material has left the grid. Although we run the simulation for 2 x 10 4 seconds, 93% of the 
stripped mass is ejected in the first 2 x 10 3 seconds, which is longer than the 1.4 x 10 3 seconds for 
the supernova to pass by the secondary. 

In Table || we compare the numerical stripped mass with an analytic estimate. The analytic 
estimate is within 10% of the numerical estimate, surprisingly close, although that changes as 
the binary separation increases. A detailed comparison of the numerical results with the analytic 
estimates is in §|. For now, we note that the analytic estimates do provide ballpark estimates of 
the stripped mass for the main sequence case. Because we use a smaller a/R ratio, consistent 
with Roche lobe overflow, our 0.15 Mq of stripped material is substantially higher than the 
0.013 — 0.052 Mq reported by Fryxell & Arnett (1981). However, we can compare our HCVc 
simulation to Fryxell & Arnett's CASE C, which has a similar a/R ratio, by scaling our stripped 
mass fraction by the new incident supernova momentum. After the scaling, our 0.022 Mq of 
stripped material increases to 0.057 Mq which is not inconsistent with the 0.052 Mq they report. 



4.4. Contamination of the Main Sequence Companion with Supernova Debris 

It is possible that some supernova material may be mixed into the envelope of the secondary 
star, either as a direct consequence of the impact or as fallback. The enhancement of heavy 
elements in the atmosphere of the secondary star may be observable 10 3 — 10 4 years (the 
Kelvin-Helmholtz timescale) after the impact or might change the course of the secondary's 
post-impact Kelvin-Helmholtz evolution by increasing the opacity of the stellar envelope. The 
question of how much supernova material is likely to be accreted by the secondary star as a direct 
result of the explosion is of particular interest in connection with the overabundance of a-elements 



observed in the atmosphere of the secondary star in Nova Scorpii 1994 (GRO J1655-40) (Israelian 



et al. 1999|) . This binary X -ray source is believed to consist of a black hole primary in close orbit 



with a subgiant secondary. Israelian et al. (1999) have interpreted the anomalous abundances in 
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the subgiant's atmosphere as contamination by fresh nucleosynthesis products ejected during a 
Type II supernova explosion that might have accompanied the formation of the black hole. 

Our calculations imply that the stellar remnant is not contaminated by the supernova material 
as a direct result of the impact. At the end of our main sequence simulation, no significant 
amount of supernova material is within the original radius of the secondary. However, a trace 
amount (1.3 x 10~ 4 M Q ) of supernova material, predominately iron-rich (> 90%) is bound by the 
secondary star, although the only supernova material in the immediate vicinity of the secondary 
star at the end of the main sequence simulation is physically located exterior to the bow shock. 
It is possible that the secondary may accrete this low-velocity iron and any additional iron mixed 
with the marginally-bound companion material when the displaced envelope, which is embedded 
within the iron layer of the supernova ejecta, begins to fall back. In addition, the slow trailing 
edge of the supernova ejecta could act as a source of iron-group elements long after the impact. 

We can make a ballpark estimate of the amount of supernova debris that can contaminate 
the secondary star by finding the amount of mass of the supernova ejecta which can not escape 
the gravitational potential of the secondary. A typical escape velocity for a 1.0 Mq main sequence 
star is ~ 600 km s" 1 . From the supernova ejecta model W7, we find that ~ 1.5 x 10 -3 M & has a 
velocity lower than the escape velocity; we can use this as a crude and conservative upper limit 
for the degree of contamination. 

4.5. Distribution in Velocity and Angle of the Material Stripped from the Main 

Sequence Companion 

One of our key results is the velocity distribution of the stripped hydrogen material, shown in 



Figure 14 a, for the HCV simulation. Also plotted in this figure is the velocity distribution of W7, 
the supernova ejecta model. The first feature to catch the eye is the clear separation in velocity 
space between the stripped material and the supernova ejecta. The velocity at the half-mass point 
in the velocity profile of the main sequence secondary is 823 km s" 1 , much less than the half- mass 
velocity of the supernova ejecta of 7836 km s _1 . Each half-mass velocity is marked with a vertical 
line and an arrow in Figure |i~4|a. The high- velocity tail shows three sharp spikes, which we identify 
with the three sharp density peaks at the composition transitions in SN W7. The composition 
groups for both the stripped material and supernova are shown; from this, it is clear that the 
bulk of the stripped material, predominately hydrogen, is embedded within the iron layer of the 
supernova ejecta. This result reinforces the belief that contaminating hydrogen and helium are 
most likely to be observed as late-time emission when the photosphere of the supernova moves to 
low- velocities, revealing the core. However, a high-velocity tail does exist, implying that a trace of 
companion hydrogen is swept up in the oxygen and silicon layers of the supernova ejecta. Given 
the upper limits on the hydrogen abundance from Type la observations near maximum light, this 
may provide a criterion for discriminating between Type la progenitor scenarios. 

In Figure |i~4| a we correct the low-velocity tail by allowing the stripped material to escape to 
infinity, free of the secondary's gravitational potential. If we assume the stripped material flows off 
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the grid at high-Mach numbers (typically M~ 5 — 30), then the specific enthalpy w = c^/(7 — 1) 
is small in comparison with the kinetic energy, and we can neglect the enthalpy in Bernouilli's 
Equation w + l/2v 2 + <f> = C, where C is a constant along a streamline. Therefore, we can estimate 



the velocity of the flow at infinity as = v — vj sc . In Figures 14 we have subtracted the escape 



velocity in quadrature. The subtraction in quadrature affects only the low-velocity tail. 



Figures 15 and [bj show the spatial distribution of the stripped mass. Figure 15 illustrates 
the distribution of the stripped material in solid angle (dM/dfi) in the HCV simulation, expressed 
in terms of the angle, a, from the downstream axis. The half-mass point of the HCV distribution 
is ~ 43°, so about half of the stripped material is confined in a cone within about ~ 43° of the 
downstream axis and the other half extends out to ~ 90°. The sudden rise along the axis for the 
HCV distribution is a consequence of the use of the differential element dil = 2ir sin ada when 



a —* 0. Figure [T(| illustrates the same spatial distribution as a polar plot. The polar angle, a, once 
again represents the angle from the downstream axis, and the radial coordinate represents the 
magnitude of dM/da. The profile clearly shows that the preferred direction is along the 43° line. 
In addition, the majority of the mass is exterior to the geometrical shadow, the shadow that the 
secondary would cast from the site of the explosion of the primary. Because of the compression of 
the "circumstellar" material along the downstream axis as the supernova sweeps around the star, 
very little stripped material is directly behind the secondary. In a 3-D calculation instabilities 
along the downstream axis are likely to mix the compressed "circumstellar" material with the 
stripped stellar material, filling in this region. 

As shown in Figure the collision of the supernova debris with the secondary will create a 
hole in the debris structure. The highest velocity debris will flow around the secondary, creating a 
cut-out region. The lower velocity material will flow around the bow shock, creating a wider hole 
in the inner ejecta. For the HCV simulation, the supernova debris will have a hole of ~ 31° in the 
outer ejecta, which will widen to ~ 40° in the inner ejecta. This corresponds to 7% and 12% of 
the surface of the supernova ejecta. Due to the supersonic nature of the flow, the hole will not 
close with time. Hence, the asymmetry in the supernova remnant may indicate the presence of a 
secondary star. 

If the stripped hydrogen were ejected spherically and uniformly mixed with the (unperturbed) 
supernova ejecta, then the mass fraction of hydrogen would be 0.6 to 0.7, almost solar, in the 
low-velocity region 10 3 km s _1 ). The hydrogen mass fraction, X(H), is shown as a function of 



supernova mass and velocity in Figures 17 and 18. The hydrogen mass fraction falls abruptly as 



the velocity increases. At 10 4 km s _1 , which is near the oxygen and silicon layers, X(H)~ 3 x 10 -3 . 
Of course, we know from the numerical simulation that the ejected material contaminates a wide 
solid angle behind the secondary. For the HCV simulation, 90% of the stripped hydrogen lies 
within 65.7° of the downstream axis which corresponds to a fractional solid angle of 29%. Figures 



17 and 18 show the expected mass fraction within the contaminated solid angle. As expected, the 



mass fraction rises dramatically. 
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4.6. The Kick Received by the Main Sequence Companion 

Table |7] lists the kicks given to the remnant companions in the main sequence and subgiant 
scenarios. We note that the red giant core, because of its extremely small size, receives a negligible 
kick. We calculate the kick from the supernova impact by finding the terminal velocity of the 
center of mass of the unstripped stellar material that remains on the grid. Although some material 
still bound to the star can flow off the grid, it is a negligible fraction (AM/M < 10~ 2 ), since the 
star reaches its terminal velocity very early in the simulation. For the main sequence simulation 
HCV, the companion receives a kick of ~ 86 km s^ 1 . Although substantial, this is still much less 
than its 227 km s _1 orbital velocity. 

In the HCV simulation the terminal velocity is reached by 5000 seconds after the explosion. 
The secondary receives a very strong kick (> 100 km s _1 ) in the first ~ 500 seconds as the 
shock passes through the center of the companion. As the shock converges and reflects in the 
back half of the companion (~ 700 seconds), the pressure rises dramatically, and the companion, 
now elongated, decelerates. By 1500 seconds, part of the shock-heated envelope has escaped the 
star and the over-pressure behind the star begins to drop. As the mass-stripping phase nears 
completion, the companion is almost spherically symmetric and close to its terminal velocity. 

Column three in Table [7| provides our ballpark estimate of the kick that we calculate by 
assuming an inelastic collision in which all the momentum of the supernova ejecta incident on 
the remnant of the companion is directly transferred to the remnant. That is, using AO^em 
and Mrem as the solid angle subtended by the remnant and the mass of the remnant, and M$n 
and Vsn as the mass and velocity of the supernova ejecta, respectively, we estimate the kick as 
v <; AQremMsnVsn/Mjiem- This is clearly an overestimate because not all of the momentum 
will be directly transferred to the remnant. The momentum transferred to the remnant will 
depend on the details of the interaction, such as the changing geometrical area as mass loss 
decreases the surface area of the companion. Our estimate is directly related to the method of 



Wheeler et al. (1975), except that we drop the logarithmic factor because it was concluded ( Fryxell 
fc Arnett 198l| ; |Taam &: Fryxell 1984j) that ablation did not substantially increase the kick. 



5. Systematic Trends with Binary Separation 

To study the effect of changing the binary separation, we performed four test cases with 
the 1.0 Mq main sequence secondary in addition to the HCV simulation: HCVa, HCVb, HCVc, 
and HCVd. The ratio of binary separation to stellar radius (a/R) is 2.57, 4.0, 6.0, and 12.0, 
respectively, with HCV having a ratio of 3.0, making it the second-closest secondary (see Table 
U). Because the incident momentum scales as the inverse square of the separation, we expect a 
significant change in the amount of stripped mass and in the kick. From Tables |6] and |7| we see 
that over this range of orbital separation the stripped mass varies more than a factor of 100, from 
0.23 Mq to 0.0018 Mq , and the kick varies by a factor of 9, from 137 km s _1 to 15 km s . 
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Figures 19 visually demonstrate the change in the strength of the impact. In the strongest 
impact, the companion experiences the strongest initial acceleration, followed by the strongest 
deceleration as the pressure rises in the back of the star. In the weakest impact, the companion 
experiences a very weak acceleration, which is not followed by a deceleration phase. In this case, 
the shock propagating through the companion is too weak to cause a pressure rise in the back of 



the star. Figures |19| illustrate the decrease in strength of the impacts with a series of 2-D images 
shown at 1000 — 1300 seconds after the supernova explosion. In the HCV simulation (Figure |l~9|a), 
the shock has already passed through the companion. The bow shock is broad, within which we see 
well-developed Kelvin-Helmholtz and Rayleigh- Taylor instabilities trailing behind the companion 
that indicate mass loss is occurring. In the HCVb simulation (Figure |l~S|b), the interaction is at 
a comparable phase, except that the shock has just converged behind the star and the trailing 
instabilities are less developed. In the HCVc simulation (Figure |i~9|c), the shock is still passing 
through the star and the bow shock is noticeably much narrower than in the stronger impacts. In 
comparison with the other cases, the impact is almost mild in the HCVd simulation (Figure |l~9|d). 
Here, the bow shock is very narrow and the shock's propagation through the star is slow. 

We use the separation study to compare our numerical stripped mass and kick with the 
analytic formulae of Wheeler et al. (1975). From Tables ||, [j], and f| we see that the stripped 
mass and kick are sensitive functions of (a/R). As shown in the tables, the analytic stripped 
mass usually overestimates the amount of mass loss by 5 — 50% for separations in the range of 
a/R = 3.0 — 4.0. For greater separations, the mass loss is severely overestimated. For smaller 
separations, the mass loss is underestimated. Our analytic estimate of the kick is calculated 
from simple conservation of momentum and does not include a logarithmic factor to represent an 
additional kick from ablation. The analytic method consistently overestimates the kick from the 
impact, by 18% — 47% in the range of our test cases. To conclude, we can say that the analytic 
estimates of the stripped mass provide a good estimate (within 10%) for binaries near Roche lobe 
overflow (a/R ~ 3), but for closer and farther separations the estimates are substantially in error. 
Furthermore, the analytic estimates of the kicks provide a ballpark value for all separations, but 
are systematically high. The kicks are below the orbital velocity of the secondary for all orbital 
separations. 

We can use the separation study to see how the average drag coefficient is likely to depend 
on the strength of the impact. The efficiency of the transfer of momentum is given by the ratio of 
stellar remnant momentum to the incident supernova momentum and can be equated to an average 
drag coefficient (Table [?]). Typical coefficients lie in the range of 0.2 — 0.5, rising systematically 
as the binary separation increases. We note that, as the binary separation increases, less mass is 
stripped from the companion and its geometrical area does not decrease as much. More of the 
incident supernova momentum can be intercepted by the remnant. In essence, the average drag 
rises for impacts in which less material is stripped. 

Despite the dramatic change in stripped mass and kick, the characteristic velocity of the 
stripped material stays between 855 km s _1 and 664 km s _1 , with the stripped material from the 
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closest binaries having the highest characteristic velocities. The half-mass velocities are: 855 km 
s _1 , 823 km s _1 , 778 km s -1 , 664 km s -1 , and 730 km s" 1 for the HCVa, HCV, HCVb, HCVc, 
and HCVd simulations respectively. The fluctuations in the high-velocity end of the velocity 
distributions seem to reflect the density variations in the initial supernova ejecta model, and like 
the characteristic velocity, decrease in strength and velocity as the binary separation increases. 

We can understand the characteristic velocity by noting that the companion's shocked 
material has a specific internal energy equal to its specific kinetic energy, which is determined by 
the speed of the shock. When the shock is strong, the envelope is heated enough that its new 
speed of sound exceeds the companion's escape velocity and the envelope is free to evaporate 
away from the star. As the shock decelerates, as it will when it climbs the steep density gradient 
of the star, less specific internal energy is deposited in the envelope and the new speed of sound 
will be lower. Eventually, the speed of sound will be lower than the escape velocity and the 
heated material will remain bound to the star. The material that escapes first will be the material 
shock-heated in the outer layers where most of the energy is deposited; this material will have the 
highest characteristic velocities. In time, the underlying layers of the envelope will evaporate off, 
leaving the marginally bound material clinging to the stellar surface. The amount of mass heated 
to sufficiently high temperatures is less for greater separations. 

The closest companions will create the widest plume of stripped material. For an explosion 
in a close binary with a hydrogen-rich companion, the contaminating hydrogen can fill a solid 
angle of 2ir steradians. A companion in a more distant binary will leave a narrower trail of 
contaminating hydrogen. Figure |2^ shows the solid angle distribution plotted as a polar plot, with 
the radius being the magnitude of dM/da and a once again being the angle from the downstream 
axis. Clearly, the closer the secondary is to the explosion, the wider in solid angle is the stripped 
material. The half-mass angles for each distribution range from 10° to 45°. As shown in Figure 



20, most of the contaminating hydrogen is exterior to the geometrical shadow. The outer extent 



of the contamination is related to the amount of mass loss, which as can be seen in Figures 19, is 
greatest for the strongest impacts. 



6. The Subgiant Companion 

Next, we describe the impact of the blast wave on a subgiant star, our HCVL scenario. As 
in the main sequence case, we neglect any changes to the structure of the star from the Roche 
potential or an optically thick wind from the white dwarf which would complicate the simulation. 
For this simulation we use a 330(p) x 565(z) cylindrical grid with the 1.13 M Q subgiant secondary 
centered on the origin. Approximately 200(p) x 100(z) zones are allocated to the secondary. We 
use the same dimensions relative to the companion's radius as the main sequence simulation: 6 
stellar radii in the p direction, 12 stellar radii in the downstream (negative z) direction and 2 
stellar radii in the upstream (positive z) direction. The same supernova ejecta model (W7) is used 
and the debris flows onto the grid in the same way. We show time frames of this simulation up to 
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2.0 x 10 4 seconds, which is longer than the companion's dynamical time of 3.5 x 10 3 seconds and 
longer than the 2.5 x 10 3 seconds that it takes for the trailing edge of the supernova to sweep past 
the back of the companion. 

For clarity, this section is subdivided into distinct topics. In §|6.1| we describe the hydrodynamic 



stages of the impact and in §6.2 we describe the stripped mass for this scenario. In §|6.3|, we discuss 



the velocity and solid angle distributions of the stripped companion material and the observational 



implications. Finally, in §16^, we discuss the kick received by the remnant of the companion. 
To avoid repetition, we focus only on the important differences between the main sequence and 
subgiant simulations. 



6.1. Hydrodynamics of the Impact on the Subgiant Companion 

The impact of the supernova W7 on the subgiant follows the same sequence of events as the 
impact on the main sequence star. The subgiant subtends almost the same solid angle as the main 
sequence star (see Table Q), so the incident momentum is almost the same in both simulations. 
Given this, it is not surprising that the character of the impact has not changed. However, because 
of the increase in physical scale, which goes as the radius of the companion star, the sequence of 
events occurs at slightly different times. Also, the compact core of the subgiant alters the progress 
of the shock through the star, and the larger envelope increases the timescale required to bring 
the remnant back into hydrostatic equilibrium. The highlights of the impact on the subgiant are 
shown in several 2-D images in Figures ^l] and f22[ 

Figure pl| a shows the impact of the supernova ejecta on the subgiant 343 seconds after the 
explosion. The leading shock front preceding the ejecta has converged along the downstream axis. 
The density jumps in the ejecta that are visible as spherical rings have already distorted the bow 
shock. More than 90% of the supernova's mass and momentum that will collide with the subgiant 
is now represented on the grid. At this early stage, only 0.06 Mq of the secondary's envelope has 
been stripped and the stripped material is still in close proximity to the secondary. 



Figure |2l| b illustrates the curvature of the shock front at 945 seconds. As with the main 
sequence companion, the shock front driven through the center of the subgiant decelerates as it 
nears the subgiant's tightly-bound core. The shock front driven into the outer envelope propagates 
so much faster that it converges on the downstream axis just after the stellar core is shocked. As 
with the main sequence companion, a reflection shock passes back through the subgiant. 



By Figure |2l|c, 2044 seconds after the explosion, the supernova shell has swept by the 
secondary. Shock-heated material can be seen streaming from the subgiant within the cavity 
shaped by the bow shock. By this time, 0.15 M & of stellar material has been stripped. The 
Rayleigh- Taylor and Kelvin-Helmholtz instabilities can be seen easily in Figure ^2[ 

At 2 x 10 4 seconds (Figure ||d), the companion is beginning to recover hydrostatic 
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equilibrium. Like the main sequence companion, the star is surrounded by an extended envelope 
of low-density material. A complicated series of shocks formed during the stellar pulsation phase 
are propagating slowly through the outer envelope. By this time, 0.17 M@ has been stripped and 
the subgiant has received a kick of 49 km s _1 from the impact. 

The pressure and density of the subgiant fluctuate in the same way as in the main sequence 
case. The central pressure first rises rapidly as the shock passes through the stellar core at 
600 — 700 seconds, then slowly falls as the star expands. As in the main sequence impact, the 
material which is now only marginally bound begins to collapse back onto the subgiant and at 
~ 4000 seconds a shock is transmitted into the outer envelope. A pulsation phase follows. 

The density and entropy profiles for the pre- and post-impact subgiant are shown in Figures 
11 and The radius of the remnant star, which we define as the radius that encloses about 90% 



of the remnant's mass, is ~ 0.9 Rq, which is a decrease of 48%. The remaining 10% of the mass 
forms a low-density extended envelope out to 5.9 x 10 6 km (~ 8 Rq). 



6.2. The Stripped Mass of the Subgiant Companion 

We determine the amount of mass stripped in the same way we did for the main sequence 
case and list it in Table |(| We include the mass stripped for two subgiant scenarios: HCVL 
and HCVLa. The HCVL binary scenario uses a 1.13 Mq subgiant which was created from a 2.1 
Mq subgiant that was artificially stripped of its outer envelope to coarsely mimic mass transfer 
prior to the supernova explosion. The HCVLa binary scenario uses the original 2.1 Mq subgiant 
with a slightly larger binary separation so that the same solid angle is subtended in order to 
ensure that both secondaries receive the same incident momentum from the blast. This scenario 
was included only as a test case. At the end of the HCVL simulation, 3 x 10 4 seconds after the 
supernova explosion (~ 5 sound crossing times), 0.17 Mq of stellar material has become unbound, 
and of this material, about 97% has left the grid by the last timestep. The mass-stripping phase 
lasts roughly 3.7 x 10 3 seconds, long enough for 96% for the mass that will be stripped by the end 
of the simulation to be stripped, and much longer than the 2.5 x 10 3 seconds required for passage 
of the supernova blast around the secondary. In comparison, the original 2.1 Mq subgiant in the 
HCVLa simulation loses 0.25 Mq. 

Table || compares the numerical stripped mass with an analytic estimate. The analytic 
calculation, which underestimates the stripped mass for both the HCVL and HCVLa simulations, 
is nonetheless within 22% of the numerical result. We ran the HCVLa simulation as a test case 
to see if the stripped mass scales directly with the companion's mass for the same subtended 
solid angle. We find that the mass loss scales almost directly with the mass of the companion. 
That is, the 2.1 Mq subgiant loses 0.25 Mq, which is 12% of its mass. In comparison, the 1.13 
M subgiant loses 0.17 Mq, which is 15% of its mass. 
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6.3. Distribution in Velocity and Angle of the Material Stripped from the Subgiant 

Companion 

Our key result for the subgiant simulation is the velocity distribution of the stripped hydrogen, 
shown in Figure |i~4| b with the SN W7 velocity distribution. The velocity profile of the material 
stripped from the subgiant shows many of the same features seen in the main sequence star 
simulation. The half-mass point for the subgiant is just slightly faster: 890 km s" 1 compared 
with 820 km s _1 for the main sequence case. Both characteristic velocities are much less than the 
half-mass velocity of the supernova ejecta of 7836 km s" 1 . Fi gure |l4] b marks the characteristic 
velocities with a vertical line and an arrow. As in the main sequence case, the bulk of the stripped 
material, predominately hydrogen and helium, is embedded within the iron layer of the supernova 
ejecta. Once again, a high-velocity tail exists which implies that a trace of companion hydrogen is 
swept up into the oxygen and silicon layers of the supernova ejecta. 



The differences between the main sequence and subgiant can be seen in Figures 14, which 



shows the asymptotic velocity profile. As discussed in §4.5, we subtract in quadrature an estimate 
for the local escape velocity that was used to distinguish between the stripped and the unstripped 
material. This approximates the velocity profile at infinity, when the stripped material is free of 
the secondary's gravitational potential. The velocity profile of the mass stripped from the subgiant 
companion shows a subtle shift to higher velocities, and the high- velocity tail shows a single peak 
near 3 x 10 3 km s _1 . The velocity profile of mass stripped from the main sequence companion 
shows three peaks, which we identify with the three composition transitions in the W7 model; the 
other two peaks in the mass stripped from the subgiant are not clearly distinguishable. 

Just as the velocity profile of the subgiant's stripped material is very similar to the velocity 
profile of the main sequence star's stripped material, the spatial distribution of the subgiant's 



stripped material shows a strong similarity to that in the main sequence case. Figure 15 
illustrates the distribution in solid angle (dM/dfi) of the stripped material in the HCVL and HCV 
simulations, expressed in terms of the angle, a, from the downstream axis. The half-mass point of 
the HCVL distribution is ~ 49°, just slightly more extended in angle than in the main sequence 



case, which has a half-mass angle of 43°. In Figure 16, the polar angle, a, once again represents 
the angle from the downstream axis and the radial coordinate represents the magnitude of dM/da. 
The material stripped from the subgiant and the main sequence companions fills almost the same 
cone, with the subgiant filling in a slightly broader one. Once again, the majority of the stripped 
mass is exterior to the geometrical shadow, the shadow that the secondary would cast from the site 
of the explosion of the primary. The collision of the supernova debris with the subgiant secondary 
will create a hole in the debris structure of ~ 32° in the high-velocity ejecta and of ~ 40° in the 
low-velocity ejecta. The asymmetry of the supernova ejecta will be slightly larger at the highest 
velocities for the subgiant model than the main sequence model. The asymmetry at low velocities 
will be almost the same for both models. As in the main sequence model, the asymmetry of the 
supernova ejecta is a direct indication of the presence of a companion star. 
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The hydrogen mass fraction is shown as a function of supernova mass and velocity in Figures 
23 and 24. The hydrogen mass fraction falls abruptly as the velocity increases. At 10 4 km 
s _1 , which is near the oxygen and silicon layers, X(H)~ 3 x 10~ 3 . Of course, we know from 
the numerical simulation that the ejected material contaminates a wide solid angle behind the 
secondary. For the HCVL simulation, 90% of the stripped hydrogen lies within 72.5° of the 
downstream axis, which corresponds to a fractional solid angle of 35%. Figures 23 and 24 show the 
expected mass fraction within the contaminated solid angle. As expected, the mass fraction rises 
dramatically. Of course, if the subgiant were of higher mass, then the hydrogen mass fraction in 
the supernova ejecta will rise accordingly. A higher mass secondary may be physically reasonable 
because our subgiant was initially 2.1 Mq before its envelope was modified to mimic mass transfer 
prior to explosion. 



6.4. The Kick Received by the Subgiant Companion 

The net kicks for the HCVL and HCVLa subgiant scenarios are listed in Table [7| along 
with the kicks for the main sequence companions. For the HCVL simulation, we see the same 
qualitative behavior as in the main sequence simulations. The secondary receives a strong kick 
in the first 700 seconds after the explosion. A brief period of deceleration occurs until ~ 1700 
seconds, and the secondary is left with a terminal velocity of 49 km s _1 by 4000 seconds. The 
numerical velocity is only 70% of the analytic value of 70 km s _1 . 



7. The Red Giant Companion 

In this section, we describe the impact of a blast wave with a red giant star, our HALGOLa 
and SYMB scenarios. These scenarios differ in the explosion model and binary separation (see 
Table ||). The SYMB scenario has a slightly larger binary separation because the mass transfer is 
by wind accretion, and its white dwarf primary explodes as a sub-Chandrasekhar Type la. The 
HALGOLa scenario has a binary separation consistent with Roche lobe overflow, and its white 
dwarf explodes as a Chandrasekhar Type la. We use the same 0.98 Mq red giant model for both 
simulations. As in the main sequence and subgiant cases, we neglect any complicating changes to 
the structure of the envelope from the Roche potential. 

We change to spherical coordinates for all the red giant simulations. The grid extends 12 
stellar radii and tt radians. We use 585(r) x 300(0) zones for the HALGOLa simulations and 
653(r) x 300(0) zones for the SYMB simulation, with 242(r) x 300(0) and 247(r) x 300(0) zones 
reserved for the red giant's envelope. The slight difference in the number of zones results from 
changing the grid to accommodate the different density and velocity profiles of the SN W7 



QNomoto et al. 1984|) and SN Hedt (|Woosley & Weaver 1994[) ejecta models. We run the SYMB 



simulation for 8.0 x 10 6 seconds and the HALGOLa simulation for 6.0 x 10 6 seconds, which is 



-33- 



only slightly more than one dynamical time, but still plenty of time for the stripped material to 
flow off the hydrodynamic grid. Because the loosely-bound envelope of the red giant is almost 
entirely stripped by the impact, there is no need to extend the calculation any further. We will 
concentrate on the SYMB simulation in this discussion, but will point out differences with the 
HALGOLa scenario when appropriate. 



For clarity, this section is subdivided into distinct topics. Section 7.1 describes the character 



of the impact and § |7.2| discusses the stripped mass for this scenario. In § |7.3| , the velocity and 
solid angle distributions of the stripped companion material and the observational implications 
are discussed. We focus only on the important differences with the main sequence and subgiant 
simulations discussed earlier. 



7.1. Hydrodynamics of the Impact on the Red Giant Companion 

Because of the lower binding energy of the envelope of the red giant (Figure ^), the impact 
of the supernova on the red giant companion is much more dramatic than the impact on the 
main sequence and subgiant companions. Unlike the main sequence and subgiant secondaries, the 
incident kinetic energy (~ 10 49 ergs) is over 1000 times the envelope binding energy (~ 10 46 ergs . 
The catastrophic effect of the blast on the weakly-bound envelope was pointed out by a number 
of authors (cf. Wheeler et al. 1975| ; |Chugai 1986| ; |Applegate fc Terman 1989| ) and demonstrated 



numerically by Livne et al. (1992). Our simulations confirm the basic result of Livne et al. (1992): 
the envelope of the red giant is stripped and the velocity of the stripped material is less than 10 3 
km s _1 . The half-mass point of our velocity distribution is 600 km s^ 1 . Like Livne et al. (1992), 
we select a binary scenario with the red giant close enough to the white dwarf to be in Roche 
lobe overflow, our HALGOLa, and a scenario with a slightly larger binary separation, our SYMB. 
Unlike Livne et al. (1992), instead of a parameterized debris model, we use supernova Type la W7 
for the blast for the HALGOLa and Type la sub-Chandrasekhar Hedt for the SYMB (Table ||). 
Figures ^5] and ^6] show highlights of the destruction of the red giant's envelope in a series of 2-D 
images from the SYMB simulation. However, Figures |2f3 focus on the propagation of the shock 



through the envelope, while Figures 26 show the entire grid 



To start the simulation, the red giant is centered at the origin, and the density and velocity 
profiles of the exploding white dwarf, scaled to 1.73 hours after the explosion, are positioned 
along the upstream axis. The impact begins at ~ 3.0 hours after the explosion. Figures p5k and 



26| a show the interaction 5.67 hours later, at 8.68 hours after the explosion. The initial impact has 
driven a shock into the companion. The shock front in the "circumstellar" medium joins smoothly 
with the shock propagating through the red giant's envelope. The supernova material has not yet 
converged on the downstream axis. 

By 14.23 hours, the bow shock has formed and the supernova ejecta have converged on the 



downstream axis right behind the star (Figures 25 b and E6p). The shock passing through the 
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companion is extremely curved. Kelvin-Helmholtz and Rayleigh- Taylor instabilities are clearly 
visible in the shock-heated material. The faint line of low-density material that is along the 
upstream axis and just outside the bow shock is a numerical artifact caused by the reflecting 
boundary. 



After 2.10 days (Figures |2q c and |2qc) , the shock moving through the envelope converges in 
the back of the star. The entire star has been shocked, and almost the entire envelope has now 
been stripped down to the degenerate core. The slight distortion of the bow shock is due to the 
numerical artifact along the upstream axis. A reflection shock will soon be created that will pass 
back through the star. 

By 5.86 days, pressure behind the bow shock has pushed it back towards the site of the 



explosion (Figure 25 d and 26 d). Although the velocity of the stripped material near the core is 
almost spherical, the bulk of the stripped material is moving downstream. In this simulation the 
bow shock is showing distortions as the density drops in the post-shock region. At this point, we 
regrid to increase the timestep and continue the simulation for an additional ~ 86.7 days, or 92.6 
days after the explosion, which is time enough for ~ 85% of the stripped material to flow off the 
grid. Only a trace amount of material, at most 4% of the original envelope (~ 0.02 Mq), is left 
clinging to the degenerate core. 

Unlike the main sequence and subgiant companions, which are done in cylindrical coordinates 
with the supernova debris added to the grid by the boundary conditions, the red giant simulations 
require the density and velocity profiles to be placed directly on the grid, off-center, as part of 
the initial conditions. The supersonic flow must propagate across radial lines and maintain its 
spherical structure. Figures ^6] show that the spherical structure of the debris is maintained. But 
does placing the spherical explosion off-center have any consequences for the amount of stripped 
mass or the velocity distribution? To test this and to address a concern that the numerical artifact 
along the axis might have unexpected consequences, we inverted the position of the supernova and 
the secondary for the HALGOLa simulation. Although there are small differences between the 
two simulations, in both cases over 0.54 Mq of the red giant's envelope is ejected by the blast. For 
the original HALGOLa simulation the half-mass point of the velocity distribution of the stripped 
material is 593 km s _1 , and for the inverted simulation the half-mass velocity is 584 km s . 

After the red giant's envelope is ejected, the degenerate core, now a single low-mass He 
pre- white dwarf, is left surrounded by a hot, entropized atmosphere. Although the blast strips 
most of the envelope, a small residual fraction (a few percent) may form an extended, hydrogen-rich 
envelope around the star. Although about ~ 75% of this envelope is within 4.4 x 10 7 km and has 
typical temperatures of 10 4 — 10 5 K, the remaining ~ 25% is cooling adiabatically, and extends out 
as far as 2.4 x 10 8 km. Typical entropies in the envelope are 35 — 40 in units of k^, 1 baryon -1 . In 
time, this extended envelope will settle onto the star, adding ~ 0.02 Mq of hydrogen-rich material 
to the hydrogen burning layer. We speculate on the fate of the low-mass He pre-white dwarf in §|9|. 
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7.2. The Stripped Mass of the Red Giant Companion 

We determine the mass lost in the same way we did for the main sequence and the subgiant 
cases. The stripped mass for the red giant simulations is shown in Table || with the results 
from the other simulations. Like Livne et al. (1992), we find that almost the entire envelope is 
removed by the impact, for all the red giant binary scenarios. For the HALGOLa simulation, the 
mass-stripping phase lasts roughly 10 5 seconds (or 1.2 days). After 6.0 x 10 6 seconds (or 69.5 
days), 0.54 Mq (98% of the stellar envelope) of stellar material has become unbound, and of this 
material, about 85% has left the grid. The stripped mass for the SYMB simulation is similar. 
The mass-stripping phase for the SYMB simulation lasts roughly 1.3 x 10 5 seconds (or 1.5 days). 
After 8.0 x 10 6 seconds (or 92.6 days), 0.53 Mq (96% of the stellar envelope) of stellar material 
has become unbound, and of this material about 85% has left the grid. 



7.3. Distribution in Velocity and Angle of the Material Stripped from the Red 

Giant Companion 



Figures 27 show the velocity distribution of the stripped envelopes; both are shown with 
the original velocity distribution of the supernova ejecta, W7 for the HALGOLa simulation and 
Hedt for the SYMB simulation. The velocity distributions from the red giant simulations, which 
have the same basic shape as the distributions from the main sequence and subgiant simulations, 
have systematically shifted to lower velocities. The half-mass velocity of HALGOLa is 593 km 
s _1 compared with 421 km s _1 for the SYMB case. 

Because the Hedt debris profile has only 78% of the momentum of the W7 profile, we expect 
the SYMB scenario to have a lower characteristic velocity. In fact, the SYMB stripped material 
has a half-mass velocity of 71% of the HALGOLa case. However, the SYMB scenario has a slightly 
greater separation. By repeating the SYMB simulation with the same binary separation as the 
HALGOLa simulation, we find that the solid angle distribution remains almost unchanged, but 
the new velocity profile is almost, but not quite, identical to the HALGOLa profile. As in the 
HALGOLa distribution, the new half-mass velocity is 593 km s . The only distinct difference 
between the velocity profiles is the slight bump in the high-velocity tail (near ~ 3.0 x 10 3 km 
s _1 in Figures |2^), which is visible in the HALGOLa case and not visible in either the old or the 
new SYMB case. We conclude from this that the slight increase in the high-velocity tail reflects 
the different supernova ejecta models used and that the shift to lower velocities in the SYMB case 
reflects the slightly greater binary separation. 

The solid angle distribution for both of the red giant scenarios is shown in Figure ^ and as a 
polar plot in Figure Although the dM/dQ profiles appear to be rising along the downstream 
axis, the amount of mass ejected in that direction is actually less than would be expected if the 
envelope had been uniformly stripped. In fact, the half-mass angles are 61° and 66°, with that for 
the SYMB scenario being slightly more extended. As in the main sequence and subgiant cases, 



-36- 



the majority of the stripped mass is exterior to the solid angle enclosed by the geometrical shadow 
of the secondary as seen from the site of the explosion. In addition, there is now a substantial 
amount of mass that has been ejected past the 90° point. 

The impact of the supernova ejecta onto the secondary creates a hole in the supernova debris 
(Figure |28| ). For the HALGOLa model, the supernova debris has a hole of ~ 34° in the outer 
ejecta, which widens to ~ 40° in the inner ejecta. This corresponds to 9% — 12% of the surface of 
the ejecta. The SYMB model has a slightly smaller hole of ~ 32° in the outer ejecta, which widens 
to ~ 37° in the inner ejecta. This corresponds to 8% — 12% of the surface of the ejecta. The 
HALGOL model clearly creates more asymmetry in the supernova ejecta at the highest velocities 
than the SYMB model or the main sequence and subgiant models. 



Figures 3C and 31 illustrate X(H), the hydrogen mass fraction in the contaminated supernova 
ejecta (assuming uniform angular mixing in the unperturbed supernova model) for the HALGOLa 
scenario. X(H) is generally higher in the red giant simulations than in the main sequence or 
subgiant simulations simply because of the larger mass ejected. However, because most of the 
hydrogen stripped from the red giant is at lower velocities, for the red giant case the hydrogen 
mass fraction at high velocities is slightly lower than in the main sequence and subgiant cases 
in that X(H) ~ 2 x 1CP 3 at velocities of 10 4 km s , compared with X(H ) ~ 3 x 10" 3 for the 
main sequence and subgiant cases. If we constrain the mixing to just the contaminated solid angle 
behind the secondary, X(H) does rise, but not dramatically, because the contaminated solid angle 
is now so much larger that this correction becomes small. For the HALGOLa simulation, 90% of 
the stripped hydrogen lies within 115.4° of the downstream axis which corresponds to a fractional 
solid angle of 71%. 



8. Upper Limit on High- Velocity Hydrogen 

Ideally, we would like a firm upper limit on the amount of hydrogen that can be mixed with 
the supernova debris at high and low velocities without being detected in the spectrum. The 



low- velocity upper limit requires a non-LTE radiative transfer study (Pinto et al. 1999). For the 
high-velocity upper limit we can use upper limits from Type la observations near maximum light. 
Although most of the stripped material is ejected at low velocities, for all the binary scenarios we 
have considered, there is a small high- velocity tail, as shown in Figure |32|. In Table §, selected 
points are tabulated. We present the velocity distribution as the total amount of hydrogen exterior 
to each velocity.^ 

Above 3.0 x 10 3 km s _1 , the left-hand edge of Figure [32], all four scenarios have less than 



4 In we discuss the effect of resolution on the high-vefocity profile. In the HCV simulation there is 15% more 
mass above 3 x 10 3 km s _1 and 26% less mass above 1.5 x 10 4 km s _1 than in the high-resolution simulation. Hence, 
at the higher velocities we conclude that changes in mass greater than ~ 30% are meaningful. 
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~ 0.014 Mq of high- velocity hydrogen. The HCVL simulation has 0.014 Mq, followed by the 
HALGOLa with 0.012 M , the HCV with 0.010 M Q , and lastly, the SYMB simulation with 
only 0.004 Mq. By 10 4 km s _1 , the amount of high- velocity hydrogen for all four scenarios has 
dropped by almost a factor of 10. The HCVL, HCV, and the HALGOLa scenarios are very similar 
up to this velocity. Because of the greater binary separation and smaller impact momentum 
and energy, the SYMB scenario has less than half the hydrogen of the other scenarios. When 
the SYMB scenario has the same binary separation as the HALGOLa scenario, it too has a 
velocity distribution like the others in this range (3 x 10 3 - 10 4 km s" 1 ). It would be difficult to 
discriminate between the main sequence, the subgiant, and the HALGOLa red giant scenarios in 
this velocity range. However, we can easily discriminate between the companions in Roche lobe 
overflow and a model in which the red giant that is too distant to be in Roche lobe overflow, our 
SYMB scenario. 

For 1.5 x 10 4 km a -1 , the HALGOLa and the SYMB scenarios have the most high-velocity 
hydrogen, the HALGOLa with 3.3 x 10" 4 M & and the SYMB with 3.0 x 10~ 4 M & . In Figure ||, 
the hydrogen distribution for the SYMB and HALGOLa scenarios converge at higher velocities, 
implying that the HALGOLa and the SYMB scenarios have almost the same amount of hydrogen 
at the highest velocities. We attribute this convergence to the momentum profile of Hedt 
(sub-Chandrasekhar model), which has more overall momentum at velocities greater than 1.5 x 10 4 
km s _1 than SN W7 (Chandrasekhar model) (Figure ||). We note that the amount of high- velocity 
hydrogen above 1.5 x 10 km s _1 could be a direct indicator of the momentum profile of the 
supernova ejecta at the highest velocities. 

Above 1.5 x 10 4 km s _1 , the HCVL and HCV scenarios have much less high-velocity hydrogen 
than in the red giant scenarios. The HCVL scenario has only 1.7 x 10 -4 Mq, which is only 53% of 
the high-velocity hydrogen of HALGOLa. The HCV has even less, only 1.3 x 10~ 4 M Q . We note 
that even if the red giant companions have more low-velocity hydrogen, they have more hydrogen 
above 1.5 x 10 km s -1 than the main sequence or subgiant models. Because the HCVL has 
more high-velocity hydrogen than the HCV scenario, in principle we can discriminate between the 
subgiant and main sequence companion at 1.5 x 10 4 km s _1 . 

The upper limit of H/Si < 2.0 x 10 -6 (relative to solar) from SN 1990M, which, assuming 
perfect mixing of the silicon in an ejecta profile like W7, corresponds to a total hydrogen mass 
of ~ 3 x 10~ 4 Mq ( pella Valle et al. 1996 ) and is indicated by a horizontal line in Figure 



Because Delia Valle et al. (1996) estimated the error in the upper limit to be a factor of 2 — 3, we 
also include horizontal lines to indicate this range. The vertical line indicates the velocity in SNIa 
W7 that corresponds to the outer 0.04 Mq of the ejecta, which they estimated to be the mass in 
the photosphere at the time of the observations. If 3 x 10 -4 Mq is a reasonable upper limit to 
the total amount of hydrogen, and if it is mixed only within the observable photosphere, then we 
can compare our mass profiles with this upper limit at the edge of the observable photosphere. 
We see from the figure that none of the binary scenarios we have considered exceeds this upper 
limit. However, if the upper limit were reduced by only a factor of 3, the red giant and subgiant 
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models would clearly exceed this limit, while the main sequence model is marginal. Because the 
supernova photosphere is non-LTE and the upper limit of Delia Valle et al. (1996) is based on 
an LTE atmosphere, the hydrogen abundance in the supernova could be much higher than their 
upper limit indicates. 

9. The Future of the Companion Star 

The post-impact evolution of the main sequence and subgiant secondaries is a separate 
calculation requiring a 2-D stellar evolution code that is capable of handling initial models that 
are slightly asymmetric in density and temperature and completely out of thermal equilibrium. 
For example, although the main sequence star is recovering hydrostatic equilibrium by the end of 
the simulation, the density in the outer envelope of the star still varies by as much as a factor of 2 
in different directions. The temperature and entropy profiles are also highly asymmetric, reflecting 
the shock history of the impact. Without attempting to follow the secondary's post-impact 
evolution in detail, we can only speculate about the secondary's future evolution. 

Immediately after the impact, the main sequence star is puffed up, much like a pre-main 
sequence star. Although the nuclear energy generation in the core has been extinguished, the 
luminosity of the residue will be dramatically brighter for this extended envelope that is out 
of thermal equilbrium. Chaboyer (1998) estimated that with the asymmetrical temperature 
distribution the luminosity could vary from 500 Lq to 5000 Lq after the impact, with a Kelvin- 
Helmholtz timescale of 1400 years to 11000 years. After thermal equilibrium is reestablished, the 
remnant will return to the main sequence along a Kelvin-Helmholtz track and then will continue 
its evolution at a rate prescribed by its new mass. Without mixing in the stellar core to refresh 
the hydrogen supply, the born-again star will appear near the middle of its main sequence lifetime. 
Using the main sequence relations R oc M 2 / 3 and L oc M 5 , the 1.0 Mq main sequence secondary 
will eventually return to the main sequence with a mass of 0.85 Mq (having lost 0.15 Mq by the 
impact), with a slightly smaller radius of 0.9 Rq, and a luminosity near 0.4 Lq. 

The subgiant puffs up after the impact and must contract to re-ignite hydrogen shell-burning, 
which will have been extinguished during the subgiant's expansion. As in the main sequence case, 
the luminosity will be extremely high during this phase. After thermal equilibrium is established 
and burning is reignited, the star will be back on a post-main sequence track with a slightly lower 
luminosity. 

After the red giant's envelope is stripped by the impact of the supernova, the degenerate core, 
now a single low-mass He pre-white dwarf, is left surrounded by a hot, entropized hydrogen/helium 
atmosphere. The structure of the core is only minimally affected by the impact because the 
converging shock was weakened by the extremely strong pressure and density gradients in the 
interior of the red giant. Although the blast strips most of the envelope, a small residual 
fraction (~ 0.02 Mq) may form an extended, hydrogen-rich envelope. Although the impact may 



-39- 



temporarily disrupt the hydrogen-burning layer, it is possible that the star may recommence 
shell burning as the envelope settles back down into a thin layer. As the degenerate core and 
hydrogen-rich envelope contracts, the star will evolve away from the red giant branch, along 
a track of constant luminosity (~ 10 3 L©) and rising effective temperature on a timescale of 
10 5 - 10 6 years for a 0.42 M© He core ( |ben fc Tutukov 1986| ; |Iben & Livio 1993| ; |Iben fc Tutukov 
1993 ). After the effective temperature rises to greater than ~ 3 x 10 4 K, the star may appear as 



an underluminous main sequence O or B star. At the end of the contraction phase, the luminosity 
begins to drop. When the luminosity drops to ~ 10 L@, it will appear as a subluminous hot star; 
an sdO star if the temperature exceeds 3.5 x 10 4 K or an sdB star if it does not ( lb en fe Livic 



1993; Green 1999). The star will then continue to cool along a standard He white dwarf cooling 



track (Benvenuto 1999). This scenario could be a possible pathway for the formation of a subset 
of single, low-mass He white dwarfs. 



10. Summary and Conclusions 

In this paper, we have presented numerous hydrodynamic simulations of the impact of a Type 
la supernova explosion with hydrogen-rich main sequence, subgiant, and red giant companions. 
The binary parameters were chosen to represent several classes of single-degenerate Type la 
progenitor models that have been suggested in the literature. All of the simulations involved 
low-mass (1.0 — 2.0 M©) companions that are close enough, or almost close enough, to be in Roche 
lobe overflow when the white dwarf primary explodes. 

We described the supernova-secondary interaction for a main sequence, subgiant, and red 
giant companions with 2-D illustrations from the simulations. The collision follows the same 
basic pattern for the main sequence and subgiant because of their similar structure. The initial 
impact of the supernova shell with the surface of the main sequence star drives a shockwave into 
the stellar envelope. A reverse shock propagates back into the ejecta. A contact discontinuity 
marks the interface between the supernova ejecta and the shocked stellar material. The shock 
propagating through the stellar envelope decelerates as it runs up the steep stellar density 
gradient. The shock front, highly curved because of the density gradients, converges in the back 
of the star. A bow shock develops in front of the companion star. After most of the supernova 
debris has passed by, the outer layers of the stellar envelope are ejected because they have been 
shock-heated to such an extent that the new speed of sound exceeds the companion's escape 
velocity. Kelvin-Helmholtz and Rayleigh- Taylor instabilities are very visible at this stage. The 
outer shocked material evaporates away from the star, embedded within the inner layer of the 
supernova ejecta. The stellar core expands and cools in response to the compression by the shock's 
passage through the core. The star has an extended, asymmetrical envelope and begins to pulsate 
as it settles back into hydrostatic equilibrium. 

The subgiant follows a similar sequence, except that its more compact core alters the progress 
of the shock through the star and its larger envelope increases the timescale required to bring 



-40- 



the remnant back into hydrostatic equilibrium. Because of its weakly-bound envelope, the impact 
on the red giant is quite different. The shock propagates through its more tenuous envelope and 
converges on its rearward side. All but a small fraction of the entire envelope is ejected, leaving a 
He pre-white dwarf surrounded by a halo of high-entropy hydrogen and helium. 

We can summarize our results by stating that as a result of the impact of the supernova 
shell the main sequence companion loses 0.15 Mq and the subgiant companion loses 0.17 Mq. 
In contrast, the red giant companions lose almost their entire envelopes, at least 0.53 Mq. The 
main sequence companion receives a kick of 86 km s _1 , the subgiant 49 km s _1 . In all cases, the 
kick received by the remnant is smaller than the original orbital velocity. Because it is too small 
to intercept more than a negligible amount of momentum, the red giant core will not receive an 
appreciable kick. 

The characteristic velocity of the stripped hydrogen (half-mass point in the distribution) is 
less than 10 3 km s _1 for all the scenarios. For the red giant case, the characteristic velocity is 
in the range 420 — 590 km s _1 , depending on the scenario. For the main sequence and subgiant 
cases, the characteristic velocities are 820 km s _1 and 890 km s _1 , respectively. With such low 
velocities, the bulk of the stripped hydrogen and helium is embedded within the low-velocity iron 
of the supernova ejecta. The region behind the stellar remnant that is contaminated with stripped 
hydrogen is always much larger than the geometrical shadow, the solid angle subtended by the 
companion as seen from the site of the explosion. For the main sequence companions the stripped 
hydrogen extends as far as 66° away from the downstream axis and for the subgiant case as far 
as 72°. The hydrogen from the red giant contaminates a much larger solid angle, as far as 115° 
from the downstream axis. These angles correspond to fractions of the sky of 29%, 35%, and 
71%. This is one distinction between the main sequence, the subgiant, and the red giant scenarios. 
The stripped hydrogen and helium is distributed over a much wider solid angle for the red giant 
scenarios. Despite the change in distribution of the stripped material, the hole in the supernova 
debris (of angular size 31° — 34°) caused by the impact with the secondary is similar for all the 
scenarios. 

If we define the efficiency of momentum transfer as the ratio of the momentum transferred 
to the stellar remnant to the incident momentum, then typical values for the main sequence 
simulations lie in the range of 0.2 — 0.5, with the highest efficiency belonging to the secondary 
with the greatest binary separation. The characteristic velocity of the stripped material stays 
within 660 — 860 km s _1 , with the material from the closest binary companions having the highest 
velocities. The closest binaries will also create the widest plume of stripped material. 

One motivation for this project is the possibility of using the velocity distribution of the 
stripped hydrogen and helium to discriminate between Type la progenitor models. However, 
to make any definitive predictions requires non-LTE radiative transfer calculations using the 
low- velocity distribution of the stripped material to determine the effect of hydrogen/helium 
contamination on the late-time supernova spectrum. Such calculations can determine the hydrogen 
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and helium line strengths and ratios of narrow emission lines that will emerge months after the 
impact, as the photosphere recedes to reveal the low- velocity supernova ejecta where the bulk of 
the stripped material is located. Unfortunately, H Q is blended with numerous Fe and Co lines, 
especially [Co III] at 6578 A, and will be difficult to identify. Because of atmospheric water 
absorption, a search for P a (1.87/im) can not easily be conducted from the ground. Searches for Pp 
(1.28/xm) can be conducted from the ground, but telluric absorption and a emission line from 
the Earth's atmosphere require careful subtraction. Also, is likely to be blended with a broad 
supernova Fe line at 1.26/j.m. Other lines to look for are He lines at 5876 A, 1.083//m, and 2.05/xm 
that will be distributed in velocity and angle like the hydrogen, but note that high-velocity helium 
is also present in sub-Chandrasekhar explosion models. With a complete study of the emergence 
of the narrow hydrogen and helium lines and the competing Fe and Co lines in the late-time 
supernova spectrum, an optimum time may be found when the hydrogen may be detected and 
used to distinguish one Type la supernova scenario from another. 

Although most of the stripped material is ejected at low velocities, all the numerical 
simulations yield a small high- velocity tail. Above 3.0 x 10 3 km s -1 , all four scenarios have no 
more than ~ 0.014 Mq of high- velocity hydrogen. The HCVL simulation has 0.014 Mq, followed 
by the HALGOLa with 0.012 M , the HCV with 0.010 M , and lastly, the SYMB simulation 
with only 0.004 Mq. Above 10 4 km s _1 , the amount of high- velocity hydrogen has dropped by a 
factor of 10 for all four scenarios. In this velocity range (3.0 x 10 3 — 10 4 km s" 1 ), we can easily 
discriminate between the companions in Roche lobe overflow and a model in which the red giant 
that is too distant to be in Roche lobe overflow, our SYMB scenario. Above 1.5 x 10 km s _1 , the 
HALGOLa and the SYMB scenario have the most high-velocity hydrogen: the HALGOLa with 
3.3 x 10~ 4 Mq and the SYMB with 3.0 x 10~ 4 Mq. We note that even if the red giant companions 
have more low-velocity hydrogen, they have substantially more hydrogen above 1.5 x 10 4 km 
s _1 than the main sequence or subgiant models. Because the HCVL has more high- velocity 
hydrogen than the HCV scenario, we can discriminate between the subgiant and main sequence 
companions at 1.5 x 10 4 km s _1 . 

The Delia Valle et al. (1996) observational upper limit from SN 1990M corresponds to a 
total hydrogen mass of ~ 3 x 10" 4 M . If this hydro gen is mixed only within the observable 
photosphere, then the main sequence, subgiant, and the red giant models are just under this upper 
limit. However, if the upper limit were reduced by only a factor of 3, the red giant and subgiant 
models would clearly exceed this limit, while the main sequence model is marginal. With this 
interpretation, the high-velocity tail of SN 1990M can not yet be used to exclude any progenitor 
models with hydrogen-rich companions. 

The impact of the supernova ejecta on the secondary star creates an oddly shaped, small gap 
in the supernova debris. The highest velocity debris flows around the secondary, creating a narrow 
gap of angle ~ 31° — 34°. The lower velocity supernova material flows through the bow shock, 
creating a wider gap in the inner ejecta of angle ~ 40°. Therefore, the hole in the supernova is 
narrow at the outside of the shell, broad at the inside of the shell, and filled with high-velocity 
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stripped material from the companion. As discussed above, most of the stripped material fills a 
wide solid angle and trails behind the supernova debris. Because the supernova flow is supersonic, 
the hole in the debris will not close with time. However, the orbital motion of the main sequence 
and subgiant secondaries will add a small component of velocity (100 — 200 km s _1 ) tangential 
to the axis of the original binary system that will slowly distort the shape of the slow-moving 
contaminated region behind most of the supernova debris. The evolution and mixing of the 
hydrogen-filled region become a 3-D problem. 

We use the angular size of the hole as an estimate of the region of asymmetry in the supernova 
ejecta caused by the impact on the companion star. The angular size of the gap in the supernova 
debris depends on the binary scenario. A typical angular size for the gap in the high-velocity 
material is 31° — 34°, which corresponds to only 7 — 8% of the ejecta's surface. The main sequence 
model has the smallest gap, the subgiant has a slightly larger gap, and the red giant has the 
largest gap. The typical angular size for the gap in the low- velocity ejecta is ~ 40°, independent of 
the model. This corresponds to 12% of the surface area of the ejecta. Because the binary scenarios 
we explore are all close enough, or almost close enough, to be in Roche lobe overflow, the degree 
of asymmetry is similar for all the models. Asymmetry in the supernova ejecta will indicate the 
presence of a companion star, but it will be difficult to use the degree of asymmetry alone to 
discriminate between main sequence, subgiant, and red giant companions. 

The asymmetry in the supernova debris could have observational consequences beyond the 
change in morphology of the supernova remnant. In the early-time spectra, subtle distortions may 
appear in the P Cygni lines that are dependent on the orientation of the hole in the supernova 
debris relative to the observer. For example, if the gap lies towards the observer, the P Cygni lines 
may have a distorted blue wing. If the gap lies at some angle away from the observer, the P Cygni 
lines may have a distorted red wing. However, because the hole in the supernova debris subtends 
such a small solid angle, its effect on the spectra will be equally small. The asymmetry may also 
affect the polarization of the supernova debris to some degree. 

The impact of the supernova shell will also have consequences for the future evolution of 
secondary star. After the impact, the main sequence star is puffed up, much like a pre-main 
sequence star. Although the energy generation rate during the impact shoots up by a factor of 
~ 1000, it immediately plummets as the star expands and cools. Hydrogen burning in the core will 
be nearly extinguished, but the luminosity will rise dramatically, to as high as ~ 5000 L Q , as the 
extended envelope relaxes back into thermal equilibrium and then follows a Kelvin-Helmholtz track 
on a timescale of ~ 10 4 yrs. The supernova is initially much more luminous with L~ 10 9 Lq. 
However, after 2 — 3 years, the supernova luminosity will have decayed enough that a very bright 
secondary might be visible. After thermal equilibrium is reestablished and the residue has settled 
back onto the main sequence, the remnant will continue its evolution at a rate prescribed by its 
new mass. The 1.0 M main sequence secondary will eventually return to the main sequence with 
a mass of 0.85 M & , having lost 0.15 M & due to the impact, with a slightly smaller radius of 0.9 
R & , and a luminosity of ~ 0.4 L & . The subgiant puffs up after the impact and must contract to 
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reignite hydrogen shell-burning, which has been extinguished during the subgiant's expansion. 
Like the main sequence case, the luminosity will be high during this phase. 

A He pre-white dwarf will be left after a supernova explosion with a red giant companion. 
Almost all the envelope will be ejected by the impact, but a residual amount of material will form 
an extended, hydrogen-rich envelope around the degenerate core. This extended envelope of ~ 0.02 
Mq of hydrogen-rich material will settle down onto the star to feed the hydrogen burning shell. 
The star will evolve away from the red giant branch at constant luminosity and increasing effective 
temperature on a timescale of 10 5 — 10 6 years. Before the luminosity drops, it may appear as an 
underluminous O or B star. When the star begins to dim and cool, it may pass through an sdO 
or sdB phase before retiring along a standard He white dwarf cooling track. This binary scenario 
could be a possible pathway for the formation of a subset of single, low-mass He white dwarfs. 

The metallicity of the secondary may increase if it accretes any supernova material as a result 
of the explosion. We do not expect that the secondary star will accrete a significant amount of 
supernova material by the direct impact of the supernova shell on the surface of the star because, 
as the shock created by the impact propagates through the stellar envelope, the shock-heated 
material escapes the secondary, carrying away any impinging supernova material. However, this 
does not necessarily mean that the secondary will not accrete any low- velocity iron-group elements 
(or oxygen and silicon if the ejecta is radially mixed) long after the initial impact. 

High resolution late-time observations focused on detecting stripped hydrogen and helium 
from the companion star are needed to verify the single-degenerate scenario. With non-LTE 
synthetic spectra, such observations may be used to discriminate between the main sequence, 
the subgiant, and the red giant companions. Asymmetry in the supernova ejecta, distorted P 
Cygni lines or line ratios, and slight polarization may also indicate an unseen companion star. 
Our theoretical study was performed in support of the quest to unravel the persistent mystery 
of the site of Type la supernova explosions. Given the emerging cosmological importance of 
such expolosions as standardizable candles and the resulting renaissance in their study, the time 
is ripe for a renewed investigation into the consequences of the variety of binary scenarios for 
thermonuclear supernovae. 
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Fig. 1. — Mass vs. radius for the main sequence, subgiant, and red giant secondaries. 

Fig. 2. — Absolute value of the exterior binding energy of the main sequence, subgiant, and red 
giant secondaries. The binding energy of each secondary is provided next to its legend. The value 
listed for the red giant secondaries is the binding energy of the envelope only. The dashed line 
indicates the binding energy profile of the red giant core. 

Fig. 3. — Frames (a)-(d) illustrate the impact of supernova Type la W7 on a 1.0 Mq main sequence 
companion at 29, 104, 179 and 329 seconds after the explosion. The color scale indicates Logiop: 
blue ~ 160.0 gm cm" 3 — ► red ~ 10~ 13 gm cm -3 . The explosion occurs off the grid at a distance 
consistent with the companion losing mass to the white dwarf primary by Roche lobe overflow. 
Frame (a) shows the companion star at the origin of the grid at 29 seconds, just prior to the 
impact of the debris. In Frame (b), at 104 seconds, the leading edge of the supernova shell has 
just collided with the surface of the companion. By 179 seconds (Frame (c)), the leading edge 
of the supernova shell converges along the downstream axis. The companion is now enveloped in 
the supernova ejecta. Frame (d), at 329 seconds, shows the shock wave that was driven into the 
companion by the impact and the bow shock just upstream of the companion. MPEG movies 



and 2-D images are posted at http: / /www.astrophysics.arizona.edu and are available via FTP from 



www.astrophysics.arizona.edu, in directory pub/marietta. 

Fig. 4. — Frames (a)-(d) illustrate the impact of supernova Type la W7 on a 1.0 Mq main sequence 
companion at 529, 755, 2029 and 2033 seconds after the explosion. The color scale is the same as 
in Figure ||. In Frame (a) the shock has just passed through the stellar core. In Frame (b) the 
shock has converged within the star along the downstream axis. Kelvin-Helmholtz and Rayleigh- 
Taylor instabilities are visible in both frames. Frame (c) illustrates the end of the mass stripping 
phase. By Frame (d), several dynamical times have passed, and the remnant is beginning to recover 
hydrostatic equilibrium. 

Fig. 5. — Frames (a)-(f) illustrate the impact of supernova Type la W7 on a 1.0 Mq main 
sequence companion in a series of cartoon images, each color indicating the dominant element 
in that region. For clarity each region is also numbered: (1) circumstellar medium, (2) stellar 
hydrogen, (3) supernova oxygen-group elements, (4) supernova silicon, and (5) supernova iron- 
group elements. The frames shown are at 104, 179, 329, 529, 755, and 2029 seconds after the 
explosion. Frames (a) and (b) show the initial impact. The oxygen, silicon, and iron layers of 
the supernova ejecta are clearly visible. Frames (c) and (d) illustrate the beginning of the mass 
stripping phase when the outermost layer of the companion's envelope is stripped and caught in the 
oxygen and silicon layers of the supernova ejecta. By Frame (e), any mass lost from the companion's 
envelope is embedded within the iron region of the supernova ejecta. Frame (f) shows the end of 
the mass stripping phase. The envelope of the companion is now extended; about 15% of the mass 
will escape the companion and follow the flow of the supernova. 



Fig. 6.— 



Momentum and density profile of SNIa W7 and SNIa Hedt. SNIa W7 is shown at 24 hrs 
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while SNIa Hedt is shown at 27.87 hrs. 

Fig. 7. — Pressure profile of the initial impact on the main sequence secondary (HCV) along the z 
axis, illustrating the propagation of the shock through the stellar core and the development of the 
bow shock. The profiles shown are at 29, 104, 179, 257, 332, 451, 601, 676, and 2030 seconds after 
the explosion and are labeled 1-9, respectively. The first and last timesteps are marked by thicker 
lines. They represent the undisturbed companion and the companion after the mass stripping is 
complete and are included for reference. 

Fig. 8. — Pressure profile of the main sequence secondary (HCV) along the z axis, illustrating 
the abrupt pressure rise in the back of the star as the stellar shock converges along the axis. As 
the shock-heated material is ejected from the star, the back pressure drops. The profiles shown 
are at 29, 751, 827, 902, 1026, 1178, 1528, and 2030 seconds after the explosion and are labeled 
1-8, respectively. The first and last timestep are marked by thicker lines. They represent the 
undisturbed companion and the companion after the mass stripping is complete. 

Fig. 9. — Pressure profile of the main sequence secondary (HCV) along the z axis as the star 
recovers hydrostatic equilibrium. The remnant is now moving down the axis with the kick received 
by the impact. The profiles shown are at 29, 2030, 5035, 10053, 15056, and 20032 seconds after the 
explosion and are labeled 1-6, respectively. The first and last timesteps are marked by thicker lines; 
they represent the undisturbed companion and the post-impact companion after several pulsations. 

Fig. 10. — Central density and temperature as a function of time for the HCV scenario (1.0 
M Q secondary) 

Fig. 11. — Pre- and Post-impact density profile of the the main sequence secondary (HCV) and the 
subgiant secondary (HCVL). 

Fig. 12. — Pre- and Post-impact entropy profile (per Boltzmann's constant per baryon) of the main 
sequence secondary (HCV) and the subgiant secondary (HCVL). 

Fig. 13. — Estimate of nuclear energy generation from the pp chains and the CNO cycle for the 
HCV scenario (1.0 M Q secondary) using a 1-D section with spherical symmetry. The compression 
and expansion of the stellar core releases ~ 10 38 ergs in the first 1000 seconds. 

Fig. 14. — A. Velocity distribution dM/dV of the stripped stellar hydrogen and helium for the 
HCV binary scenario at the end of the simulation, along with the velocity distribution of SNIa W7. 
The velocity at the half-mass point of the HCV distribution, as indicated by the vertical line and 
the arrow, is 823 km s _1 . In contrast, the velocity at the half-mass point of the SN W7 distribution, 
also indicated in the plot, is 7836 km s _1 . B. Velocity distribution dM/dV of the stripped stellar 
hydrogen and helium for the HCVL binary scenario at the end of the simulation, along with the 
velocity distribution of SNIa W7. The velocity at the half-mass point of the HCVL distribution, 
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as indicated by the vertical line and the arrow, is 887 km s . 

Fig. 15. — The solid angle distribution of the stripped envelope of the main sequence and subgiant 
simulations HCV and HCVL. The half-mass angle of each distribution is listed. The dashed lines 
indicate the distribution of supernova material for both simulations. The hole in the supernova 
debris created by the impact is visible between 0° — 40°. 

Fig. 16. — The magnitude of the angular distribution (dM/da) is plotted as the radial coordinate 
for the HCV and HCVL binary scenarios. The angular coordinate is the angle from the downstream 
axis. The geometrical shadow is shown for both simulations. 

Fig. 17. — Hydrogen mass fraction of the contaminated supernova ejecta for the HCV scenario 
(1.0 Mq main sequence secondary) vs. the mass of the original W7 supernova model. The solid line 
indicates the hydrogen mass fraction assuming that the stripped material is uniformly mixed in 
angle with Type la model W7. The mass fractions of the O, Si, and Fe-group elements are shown 
for comparison. Because the stripped material in this simulation contaminates only the supernova 
ejecta within 65.7° of the downstream axis, the dashed line indicates the hydrogen mass fraction 
assuming mixing only in this region. 

Fig. 18. — Hydrogen mass fraction of the contaminated supernova ejecta for the HCV scenario 



(1.0 Mq main sequence secondary) vs. velocity. See the caption for Figure 17 



Fig. 19. — Frames (a)-(d) illustrate the impact of supernova Type la W7 on a 1.0 Mq main 
sequence companion for different binary separations: simulations HCV, HCVb, HCVc, and HCVd. 
The color scale indicates Logiop: blue ~ 160.0 gm cm -3 — > red ~ 10 -13 gm cm -3 . This color scale 
is the same as in Figure ^. Frame (a) shows the impact for a ratio a/R = 3.0 (binary separation 
over stellar radius) at 1031 seconds after the explosion. This is our HCV simulation. Frame (b) 
shows the impact with a/R = 4.0 at 1336 seconds after explosion, the HCVb test case. Frame (c) 
shows the impact with a/R = 6.0 at 1197 seconds, and Frame (d) with a/R = 12.0 at 1308 seconds. 
The timesteps are increasing primarily because the length of time needed for the supernova ejecta 
to reach the secondary is increasing. 

Fig. 20. — The magnitude of the angular distribution (dM/da) is plotted as the radial coordinate 
for the main sequence simulation HCV and test cases HCVa, HCVb, and HCVc, which differ 
only in their distance from the supernova explosion. The half-mass angles of the distributions are 
45.1°, 42.7°, 35.8°, 24.7°, and 9.9° for the HCVa, HCV, HCVb, HCVc, and HCVd simulations, 
respectively. 

Fig. 21. — Frames (a)-(d) illustrate the impact of supernova Type la W7 on a 1.13 M & subgiant 
companion at 343, 945, 2044 and 20044 seconds after the explosion. The color scale indicates 
Logiop: blue ~ 2.7 x 10 3 gm cm -3 — » red ~ 10~ 13 gm cm" 3 . The explosion occurs off the grid at 
a distance consistent with the companion losing mass to the white dwarf primary by Roche lobe 
overflow. Frame (a) shows the subgiant and the incoming supernova shell 343 seconds after the 
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explosion. Only 0.06 Mq of the companion's envelope has been stripped by this time. In Frame (b) 
at 945 seconds, the shock has just propagated through the subgiant's core and is about to converge 
in the back of the subgiant. Frame (c) shows the subgiant in the last stages of mass loss. The 
subgiant has lost 0.15 Mq of its envelope which can be seen streaming away from the subgiant. By 
Frame (d), 20044 seconds after the explosion, the remnant of the subgiant is recovering hydrostatic 
equilibrium. 

Fig. 22. — Impact of supernova Type la W7 on a 1.13 Mq subgiant companion at 1958 seconds 
near the end of the mass stripping phase. The Kelvin-Helmholtz instabilities, which help the hot 
material to escape the surface, are easily seen in this image. The color scale is the same as in Figure 



21 



Fig. 23. — Hydrogen mass fraction of the contaminated supernova ejecta for the HCVL scenario (1.1 
Mq subgiant secondary) vs. the mass of the original W7 supernova model. The solid line indicates 
the hydrogen mass fraction assuming that the stripped material is uniformly mixed in angle with 
W7. The mass fractions of the O, Si, and Fe-group elements are shown for comparison. Because 
the stripped material in this simulation contaminates the supernova ejecta within only 72.5° of the 
downstream axis, the dashed line indicates the hydrogen mass fraction assuming mixing only in 
this region. 

Fig. 24. — Hydrogen mass fraction of the contaminated supernova ejecta for the HCVL scenario 



(1.1 M & subgiant secondary) vs. velocity. See the caption for Figure |23. 



Fig. 25. — Frames (a)-(d) illustrate the impact of supernova Type la Hedt on a 0.98 Mq red giant 
companion at 8.7 hrs, 14.3 hrs, 2.1 days and 5.9 days after the explosion. Only the region around 



the red giant is shown to emphasize the interaction. The full images are shown in Figures 26 



Fig. 26. — Frames (a)-(d) illustrate the impact of supernova Type la Hedt on a 0.98 Mq red 
giant companion at 8.7 hrs, 14.3 hrs, 2.1 days and 5.9 days after the explosion. In these images 
the supernova can been seen exploding on the grid. These are the full images corresponding to the 
images in Figures |25|. 

Fig. 27. — A. Velocity distribution dM/dV of the stripped stellar hydrogen and helium for the 
HALGOLa binary scenario at the end of the simulation, along with the velocity distribution of 
SNIa W7. The velocity at the half-mass point of the HALGOLa distribution, as indicated by the 
vertical line and the arrow, is 593 km s . In contrast, the velocity at the half-mass point of the 
SN W7 distribution, also indicated in the plot, is 7836 km s . B. Velocity distribution dM/dV 
of the stripped stellar hydrogen and helium for the SYMB binary scenario. The velocity at the 
half-mass point of the SYMB distribution, as indicated by the vertical line and the arrow, is 421 
km s _1 . The velocity at the half-mass point of the SN Hedt distribution, also indicated in the plot, 
is 10000 km s _1 . 



Fig. 28.— 



Solid angle distribution of the stripped envelope of the red giant simulations HALGOLa 
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and SYMB. The dashed lines indicate the distribution of supernova material. The gap between 
0° — 40° illustrates the hole in the supernova ejecta created by the impact. 

Fig. 29. — The magnitude of the angular distribution (dM/da) is plotted as the radial coordinate 
for the red giant simulations (HALGOLa and SYMB). 

Fig. 30. — Hydrogen mass fraction of the contaminated supernova ejecta for the HALGOLa 
scenario (0.98 M red giant secondary) vs. the mass of the original W7 supernova model. The 
solid line indicates the hydrogen mass fraction assuming that the stripped material is uniformly 
mixed in angle with Type la model W7. The mass fractions of the O, Si, and Fe-group elements 
are shown for comparison. Because the stripped material in this simulation contaminates the 
supernova ejecta within only 115.4° of the downstream axis, the dashed line indicates the hydrogen 
mass fraction assuming mixing only in this region. 

Fig. 31. — Hydrogen mass fraction of the contaminated supernova ejecta for the HALGOLa 
scenario (0.98 M & subgiant secondary) vs. velocity. See the caption for Figure |30[ 

Fig. 32. — Mass of stripped hydrogen with velocity > V for the HCV (1.0 main sequence), HCVL 
(1.13 Mq subgiant), SYMB and HALGOLa (1.0 Mq red giant) simulations. To compare the mass 
profiles to Type la supernova observations near maximum light, we indicate with a horizontal line 
the Delia Valle et al. (1996) upper limit of 3 x 10~ 4 M of hydro gen in the spectrum of SN 1990M. 
This assumes perfect radial mixing of the silicon and hydrogen. Because Delia Valle et al. (1996) 
estimated the error in the upper limit to be a factor of 2 — 3, we also include horizontal lines to 
indicate this range. Delia Valle et al. (1996) estimated that the size of the photosphere was 0.04 
Mq at the time of the observation. The velocity corresponding to this mass is estimated from SNIa 
W7 and shown in the figure as a vertical line. 
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Tablc 1. Single Degenerate Type la Models 



Single Degenerate SNIa 


Secondary 


Mass Transfer 


Hydrogen Cataclysmic Variables 


Main Sequence/Subgiant 


RLOF, H 


Hydrogen Cataclysmic-Like Variables 


Main Sequence/Subgiant 


RLOF, H 


Symbiotic Stars 


Red Giant 


Wind, H 


Hydrogen Algols 


Red Giant 


RLOF, H 


Helium Cataclysmic Variables 


He Star 


RLOF, He 


Helium Algols 


He Giant 


RLOF, He 



Table 2. Secondary Models 



Evolutionary Stage 


M (M ) 


R (km) 


tdyn (s) a 


tcross (s) b 


Main Sequence 


1.017 


6.80 x 10 5 


1.5 x 10 3 


2.6 x 10 3 


Subgiant 


2.118 


1.62 x 10 6 


4.0 x 10 3 


5.8 x 10 3 


Subgiant, mass loss 


1.132 


1.22 x 10 6 


3.5 x 10 3 


5.5 x 10 3 


Red Giant c 


0.977 


1.19 x 10 8 


6.4 x 10 6 


5.5 x 10 6 



a Dynamical time, tdyn = \ (G(p)) 1 '" 2 

Sound crossing time, t cross — R/(c a ), integrated over star 
c Softened <j>(r) = -GM c /(r + r c ), M c = 0.424 M , R c = 5 x 10 5 km 
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Table 3. Supernova Ejecta Models 



Model 


M (M ) 


E (10 51 ergs) 


P (M V c ) a 


(V) h 




SNIa W7 d 


1.378 


1.233 


1.175 


8.527 x 10 3 


7.836 x 10 3 


SNIa Hedtbll c 


0.899 


1.079 


0.912 


1.014 x 10 4 


1.006 x 10 4 



a Momentum of supernova ejecta, specified in units of Mq and divided by V c = 10 9 cm s 

b Average speed of the supernova ejecta, km s _1 

c Speed of the supernova ejecta at the half-mass point, km s _1 

d |Nomoto et al. 1984 

e [Woosley fc Weaver 1994 



Table 4. Simulations 



Simulations 




M (M ) 


a (km) 


a/R 


SN Model 


HCV 




1.017 (MS) 


2.04 x 10 6 


3.00 


W7 




Additional Separations 










HCVa 


1.017 (MS) 


1.75 x 10 6 


2.57 


W7 




HCVb 


1.017 (MS) 


2.72 x 10 6 


4.00 


W7 




HCVc 


1.017 (MS) 


4.08 x 10 6 


6.00 


W7 




HCVd 


1.017 (MS) 


8.16 x 10 6 


12.00 


W7 


HCVL 




1.132 (SG) 


3.39 x 10 6 


2.78 


W7 




No Mass Loss 










HCVLa 


2.118 (SG) 


4.50 x 10 6 


2.78 


W7 


HALGOLa 




0.977 (Red Giant) 


3.00 x 10 8 


2.52 


W7 


SYMB 




0.977 (Red Giant) 


3.76 x 10 8 


3.16 


Hedt 



Table 5. Pre- and Post-Impact Central Values for HCV Companion 



Initial Final Fractional Change 



Pressure (dynes cm" 2 ) 2.34 x 10 17 9.41 x 10 16 0.60 

Density (gm cm" 3 ) 1.54 x 10 2 7.95 x 10 1 0.48 

Temperature (K) 1.56 x 10 7 9.68 x 10 6 0.38 

Entropy (k^ 1 baryon" 1 ) 1.10 x 10 1 1.35 x 10 1 0.23 
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Table 6. Stripped Mass 



Scenario 


Time (s) 


AM (M ) a 


AM (M ) b 


HCV 


2.0 x 10 4 


0.15 


0.16 


HCVa 


2.0 x 10 4 


0.23 


0.20 


HCVb 


2.5 x 10 4 


0.074 


0.11 


HCVc 


3.0 x 10 4 


0.022 


0.063 


HCVd 


3.0 x 10 4 


0.0018 


0.021 


HCVL 


3.0 x 10 4 


0.17 


0.14 


HCVLa 


3.0 x 10 4 


0.25 


0.21 


HALGOLa 


6.0 x 10 6 


0.54 (98%) 


0.55 (100%) 


SYMB 


8.0 x 10 6 


0.53 (96%) 


0.55 (100%) 



a Numerical Calculation of Stripped Mass 

b Analytic Estimate of Stripped Mass, method of Wheeler et al. 1975 



Table 7. Kicks 



Scenario 


M rem (M ) a 


V b 

v rem 


V c 


Verb d 


P /P e 

x rem/ - 1 in 


HCV 


0.867 


85.7 


109.8 


227.1 


0.228 


HCVa 


0.787 


99.2 


137.1 


245.2 


0.176 


HCVb 


0.943 


61.5 


73.2 


196.7 


0.316 


HCVc 


0.995 


35.2 


41.7 


160.6 


0.429 


HCVd 


1.015 


10.1 


14.8 


113.5 


0.503 


HCVL 


0.962 


49.4 


70.4 


172.1 


0.125 


HCVLa 


1.868 


54.4 


65.4 


126.5 


0.267 



a Numerical Calculation of the Mass of the Stellar Remnant 
b Numerical Calculation of the Kick of the Stellar Remnant, km s _1 
c Ballpark Estimate of Kick from Inelastic Collision, km s _1 
d Orbital Velocity, km s _1 

c Ratio of Stellar Remnant Momentum (P re m) to Incident Supernova Momentum (Pi„) 
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Table 8. Mass of Stripped Hydrogen with Velocity > V at four chosen points: 
Va = 3.0 x 10 3 km s -1 , Y B = 1.0 x 10 4 km s -1 , Y c = 1-2 x 10 4 km s -1 , and 

V D = 1.5 x 10 4 km s _1 . 



Scenario 






A 










B 










C 










D 






HCV 


9.7 


X 


io- 


-3 


M 


1.3 


X 


10" 


-3 


M e 


6.9 


X 


10' 


-4 


M 


1.3 


X 


io- 


-4 


M 


HCVL 


1.4 


X 


10" 


-2 


M 


1.4 


X 


10" 


-3 


M Q 


7.8 


X 


10' 


-4 


M 


1.7 


X 


10" 


-4 


M 


HALGOLa 


1.2 


X 


10" 


-2 


M 


1.4 


X 


10" 


-3 


Mq 


8.9 


X 


io- 


-4 


M 


3.3 


X 


10" 


-4 


M 


SYMB 


4.3 


X 


10" 


-3 


M 


6.2 


X 


10" 


-4 


Mq 


4.8 


X 


io- 


-4 


M 


3.0 


X 


10" 


-4 


M & 
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